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Abstract 


A detailed investigation of the sooting structure of an overventilated 
laminar axisymmetric ethene/air diffusion flame under atmospheric pressure is 
presented. A thermophoretic sampling technique using a variety of probes is 
employed to obtain data on the morphological character of soot particles 
collected from various locations of the flame. The morphological features of 
the particles examined by Transmission Electron Microscopy provide not only 
qualitative information on _ particle agglomeration, surface growth, and 
Oxidation but also quantitative data on primary particle size as a function of 
flame coordinates. An increasing primary particle size with height above the 
burner mouth proved the dominance of surface growth reactions on the 
morphology of particles transported in the lower and intermediate portions of 
the flame. The surface growth reaction rate throughout these regions was 
calculated from a series of observations without any assumption on the form 
of the particle size distribution. The primary particle size was found to 
decrease as the particles move to the upper regions of the flame probably 
because of oxidation. A theoretical model examining the aerosol dynamics 
when particle formation, coagulation, and surface growth are present is also 
discussed. The model uses a time varying bimodal representation of the 
aerosol size distribution function with a prescribed logarithmic-normal form. 
The model can be used to study the physical/chemical processes acting upon 
the particulate phase in various aerosol systems in which particle generation 
mechanisms are present. The thermophoretic method, combined with existing 
data obtained by others through laser-light scattering and _ extinction 
measurements on the same flame, is used to describe the physical aspects of 
soot formation and transport. The combined observations of the two methods 
led to the conclusion that the soot particles are formed locally low in the 
flame in regions of high temperature near the reaction zone. The two 
experimental methods together with the theoretical model form a powerful unit 
used to gain a better understanding of the soot formation processes related to 
a diffusion flame environment. 
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CHAPTER 1 


INTRODUCTION 


Soot emission from various combustion appliances is a phenomenon 
observed very frequently. The light of a candle flame has its origin to the 
presence of incandescent soot particles. The"efficiency’ of* an ~ internal 
combustion engine, particularly of a diesel motor, depends on the suppression 
of the soot formation processes. It is therefore desirable to control the soot 
formation mechanisms in order to avoid poor combustion and a loss of 
efficiency in practical combustion systems. However the presence of soot is 
desirable in some technical processes for a more efficient heat transfer, since 
the ability of the soot particles to radiate promotes the heat transport. In 
these systems subsequent oxidation of the particles is necessary in order to 
avoid pollution. The production of carbon black from the pyrolysis of 


hydrocarbons is another process where a maximum yield of soot is desirable. 


The above discussion explains the intense research activity related to 
the soot formation processes during the last twenty years, mainly focused on 
the study of laboratory scale systems. These studies produce valuable 
-information on the soot formation mechanisms in simple combustion systems 
and can be used to study the more complicated phenomena related to large 


scale fires. 


The conditions under which soot is formed vary, and depend on 
numerous factors. Generally the generation of soot in combustion systems 
occurs in a temperature range between 1500 and 2500 K. The conversion of 
carbon available in the fuel consumed to soot, varies for different 


systems. For a laminar axisymmetric ethene/air diffusion flame under 


atmospheric pressure the conversion®® is 10-15%. However at elevated 


pressures the soot yield can be appreciably higher (~50% for the laminar 
diffusion flames studied by Flower and Bowman?®?’), The soot particles look 
black and consist mainly of carbon and up to 10% per mole hydrogen. Newly 
formed particles in flames have a larger hydrogen content and display a 
considerably higher reactivity when compared to the old particles*®. The 
macroscopic density of old soot has been found?%*46? to be in the interval 1.8 
- 2 gr/cm*. The density of incipient soot may be substantially lower because of 


the high hydrogen content. 


The soot particles inside a combustion environment form an aerosol 
which can be characterized by the following quantities: the volume 
concentration of the particulate phase 9, called particle volume fraction 
(cm?/cm?), the particle number concentration N (cm7%), and the particle size D 
(nm). Various studies concentrate on the determination of those parameters 


which uniquely define the dispersed phase. 


In the various combustion systems the time available for the 
development of the soot formation processes is of the order of a few 
milliseconds. During this period a certain sequence of events takes place in a 
specific order. The first stage, called particle inception or particle formation, 
incorporates the generation of soot "“nuclei" from the fuel molecules via 
oxidation and pyrolysis products (molecules, radicals, ions). Chemical 
kinetically controlled reactions initially create the precursor species needed to 
form the first soot particles, called incipient particles. Various theories 
about the nature of the soot precursors in flames have been developed so far. 


The most likely precursors are considered to be the unsaturated hydrocarbons, 
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like acetylene and its higher analogs (C, H,), as well as polycyclic aromatic 
hydrocarbons. Other theories suggest that flame generated ions or radicals 
are also major species involved in the soot formation processes. The precursor 
species are then combined through a mechanism that is not well understood to 
form the incipient soot particles. These particles are small*#® (1~3 nm), and 
are created in large numbers thus forming the soot aerosol hosted by the 
combustion gases. Even though the first soot particles are created locally in 
high number concentrations they are associated with a very low soot loading 


because of their small size. 


The second stage involves the growth of the existing particles 
_(increasing D) through two different mechanisms. Particle coagulation (when 
particles collide and coalesce) which is representative of the early stage of soot 
formation, or particle agglomeration (when chain-like particles are formed) 
representative of the stages thereafter. Both coagulation and agglomeration 
result in a decrease of the number concentration N while the volume fraction 
@ remains constant. The other mechanism, called surface growth, involves the 
gas phase mass deposition on the surface of the existing particles through 
heterogeneous reactions. Surface growth results in an increase of $ while N 
remains unchanged. The mechanisms of growth are mainly responsible for 
the bulk of soot generated during combustion. Agglomeration and surface 
growth may coexist and their relative contribution to the formation of the 
macroscopic fluffy soot flocs visible in the atmosphere varies accordingly to 


the conditions present in the combustion environment. 


The stages of soot formation and growth finally terminate when the 


oxidation mechanisms dominate. The third stage thus involves the burning of 
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the particles in the presence of oxidizing species to form gaseous products such 
as. CO, CO)ycetc: The intensity of the oxidative attack determines the 
eventual emission of soot particles from the combustion device. This emission 
highly depends on the temperature** and velocity fields inside the system as 


well as the presence of additives*®»®5 


which can alter its behavior drastically. 
Most soot studies so far have concentrated on simple combustion systems and 
in particular laboratory scale flames. In the premixed flames both fuel and 
oxidant are mixed prior to the combustion process. Diffusion flames are those 
in which the mixing of fuel and oxidant occurs in the combustion region 


through molecular or turbulent diffusion processes. The premixed flames have 


been thoroughly24:6,7,38-43,47,48,76,98 


investigated, and the soot formation and 
depletion mechanisms are better understood at this instant. The mechanisms of 
soot formation and growth may be quite different from system to system, and 
most of the details of these mechanisms are still not well established. This is 
clearly demonstrated by a series of recent review  articles!,45.53,55,82,96 
summarizing the results of numerous studies on soot formation. The intense 
interest of the scientific community for a better understanding of the soot 


formation processes is illustrated by the long list of 388 references in one* of 


the above mentioned reviews. 


_ 2 iesisoCiective 


The first attempts to understand the reninnie.. dP diffusion flames 
were made by Burke and Schumann? who developed a theory to predict the 
ideal flame front, based on a model of equal interdiffusion of both fuel and 
oxidizer. Their simple theory predicted the heights of a variety of 


axisymmetric flames very successfully, using the appropriate values of 
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diffusion coefficients. However the visible flame height is not always the 
same, depending on the overall production of soot as well as its depletion by 
oxidation. Most experimental work since then has concentrated on the 
investigation of soot formation, growth and depletion mechanisms in both 
premixed and diffusion flames. Although numerous studies have’ been 
completed so far, many details related to soot production and oxidation in 


flames are not yet well established. 


‘The present study attempts a detailed investigation of the sooting 
structure of an overventilated laminar axisymmetric ethene/air diffusion flame 
under atmospheric pressure conditions. A wealth of existing data provides us 
with valuable information on the temperature and velocity fields of an 


5 


identical flame®® analyzed at the Centers of Fire Research and Chemical 


Engineering of the National Bureau of Standards (NBS). A series of optical 


tests®2-86 


performed there, produced information on the structure of the soot 
field throughout the coannular ethene diffusion flame. However these tests 
proved inadequate to provide information on the size distribution of the soot 
particles, and to assess the relative contributions of various mechanisms 


(particle inception, agglomeration, surface growth, oxidation) affecting the soot 


particle morphology. 


A thermophoretic sampling technique presented in this work provides 
data on the morphological character of soot particles corresponding to various 
locations of the flame. Additional information on the contribution of surface 
growth is also obtained through the analysis of soot samples collected 
thermophoretically from the flame environment. A theoretical model of soot 


aerosol dynamics is also developed and used to interpret the experimental data. 
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The two experimental methods together with the theoretical model form a 
powerful unit, used to gain a better understanding of the soot formation 


processes related to a diffusion flame environment. 


The current work begins with a general discussion on soot particle 
diagnostic methods and a brief review of the data obtained at NBS through 
the laser light scattering and extinction measurements. A series of boundary 
layer studies on soot particle transport related to the thermophoretic method is 
presented in the next section, followed by a detailed description of the 
experimental apparatus. The experimental results are then presented with an 
extensive evaluation of the data and a critical discussion. The theorétical 
model is developed in the following section and is subsequently used to 
evaluate the data obtained experimentally. The last section presents the 
conclusions of the study and recommendations for future work. The treatise 


is concluded with the appendices and other related material. 


CHAPTER 2 


SOOT FORMATION IN A LAMINAR AXISYMMETRIC 
ETHENE/AIR_ DIFFUSION FLAME. PREVIOUS STUDIES 


The interest in the formation, evolution and partial or complete 
oxidation of soot in flames resulted in the development of a wide variety of 
sampling methods. The two most widely used are the "in situ" methods 
examining the particles in location, and the direct sampling methods which 
remove the particles from their environment. The scattering of a light beam 
has been used initially to analyze the soot field in a flame environment**”!, 
The more accurate laser light scattering and transmission tests were 


d5.14-17,44,49,51,52 for an "in situ" analysis of the soot field 


subsequently employe 
in flames. These methods have considerable appeal because they are 
nonintrusive and allow observation of the soot field without intervening in the 
chemical and physical processes which are present. However a weak point of 
these methods is their inability to distinguish between single and agglomerated 
particles that cause the scattering of the light. Direct sampling of the particles, 
an intrusive technique, has appeared’ 'in®*several’* forms so * far’ but**it ‘has 
received more limited attention because of questions related to the validity of 
the morphological observations. Some of the direct sampling methods are based 
on the fact that soot deposition rates on cold targets immersed in flames are 
dominated by particle thermophoresis, i.e. particle drift down the temperature 
gradient. According to this principle, soot particles have often been collected 
on cold surfaces*?4:76,80,88,95 for subsequent analysis by electron microscopy. 


12,71,78 to collect the 


Another method consists of using a sampling orifice 
particles from the interior of the flame and subsequently analyze them through 


electron microscopy techniques. This method however involves a few steps 


which aim at the separation of the sampled particles from the combustion 
gases withdrawn from the flame. These steps may affect the particle 
morphology thus causing a degree of uncertainty of the results obtained using 
this method. The molecular beam-sampling method has also been used 
extensively?4#%% but its use is limited only to the investigation of low pressure 
flat flames. None of the above sampling methods has clearly an edge over 
the others. A possible combination of these techniques has a greater potential 
to produce more information related to the soot formation processes in a flame 


environment. 


2.1 The ethene laminar diffusion flame, Optical studies 


The combustion system studied in the current work is a coannular 
laminar ethene/air diffusion flame under atmospheric pressure conditions. The 
burner consists of two concentric brass tubes of 11.1 mm and 101.6 mm i.d. 
The fuel flows through the inner tube while the ventilating air flow passes 
through the annulus between the two tubes. A detailed description of the 


burner is given in Section 3.4. 


The flow rates adopted throughout this work for the ethene and the 
air are 3.85 and 713.3 cm/s respectively. These flow rates result in a 
laminar flame of definite shape and height (88 mm) shown in Figure 2.1. It is 
Clear that the flame has an axisymmetric structure characterized by a 
cylindrical coordinate system (r,z). The; lowest spart ‘tofiisthe wflames\ is 
surrounded by a blue zone which extends up to a height of 3-5 mm. The rest 
of the flame appears as an envelope characterized by the familiar soot 


luminosity and a strong absorptivity. The luminous flame envelope closes at 
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the top which is a general characteristic of flames that release no soot to the 


environment (non smoking flames). 


The structure of the soot field of the coannular ethene diffusion 


flame has been previously investigated using laser light-scattering and 


83-86 


extinction measurements The apparatus used is described in ref. 86. A 


4W argon ion laser operating at the 514.5 nm laser line was used as the light 


source for the scattering and extinction measurements. The laser was 
Operated at the 488.0 nm line for the fluorescence experiments. The beam 
diameter was approximately 0.2 mm with a radial control of 0.127 mm. In 
all the experiments the incident beam was vertically polarized. We should 


mention at these point that the optical measurements were performed using a 
metal chimney instead of a plastic one that was used throughout our 
experiments. This change though is estimated to have only a slight effect on 
the structure of the flame, mainly because of the distance of the chimney 


walls from the flame (~8 burner radii). 


For particles@in: thes Ravylershwregime +(D<4/10° D Sparticle’diameter} 
wavelength of incident beam), the value of the volumetric scattering cross 


section for incident and scattered rays polarized in the vertical plane is given 


by x 
F(m) rays 
Q, = = Ee) N P(D) D® dD (2.1) 
O 
where 
2 
ts flak tony 
F(m) = ; 
m*+ 2 


~ 


and P(D) is the particle size probability function [ f P( Died D: p= 1), m the 
O 

complex refractive index of the particulate material and N the total number 

of particles per unit volume. The product P(D) dD expresses the fraction of 


all particles contained in the interval D to D+dD. The extinction coefficient 


in the Rayleigh limit is given by 


1k 
Kou = 5 Ef) N P(D) D® dD (2.2) 
ie) 
where 
m?-1 
E( ih) = -Im 
m? + 2 


A parameter thoroughly used in the optical method, called generalized 


mean diameter or moment ratio, is defined by 


| P(D)DP dD 
O 


| P(D)D4 dD 
0 


The particle volume fraction $¢, defined as the volume concentration of the 


particulate material inside the host gas, is related to K,, through 
des g2-5t sg 3 (2.4) 
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From Eqs. (2.1) and (2.2) using the definition of Eq. (2.3) we obtain 


Ly. 
4? Seay MQ,7 we 
Dig weed wr - (2.5) 
F(m) KY 
If 
3 
f vss (2.6) 
No ‘ 
Ds 
J Tl 3 
and since Ps= ap ND,p 
we obtain the following expression for the number concentration 
69 
N = ——. fy Gah 
11D 55 


All three Eee (2.4), (2.5) and (2.7) are valid independently of the 
form of P(D), provided that the particles belong to the Rayleigh regime. An 
assumption for the value of m is required in order to calculate $ and Daa 
from Eqs. (2.4) and (2.5) respectively. An additional assumption is required 
for the value of fy, in order to obtain the values of N from Eq. (2.7). The 
assumed value of fy is directly related to the form of the distribution as we 
can see from Eas. (2.3) and. (2.6). The optical measurements discussed here 
used a value of m = 1.57 - 0.56i. A different value of m would affect the 


values of N by a significant factor, while the values of $? and D,, would be 
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less influenced®*®>. The value of fy, was taken to be - equal” toe e 
corresponding to the self-preserving distribution function®’, which is the 
asymptotic form of the particle distribution function for an aerosol undergoing 
free molecular coagulation by Brownian motion in the absence of particle 


formation or surface growth. 


Soot volume fraction, particle size and number’ concentration 
measurements for the ethene flame are shown in Figure 2.2 for two heights 
(z=15 and 50 mm). Only the labeled profiles are of interest for our work, 
since the rest correspond to other flames studied by Santoro et al®®, The 
soot volume fraction at each height achieves a maximum value at a radial 
location r, defining the annular region. This maximum value steadily 
increases up to a height of 40 mm, reaches a maximum between z=40 mm and 
z=50 mm and then decreases probably through oxidation. The central region 
of the flame (r=0) is first observed to contain particles at an axial location 
which is significantly higher than that observed for the annular region. The 
regions of maximum @ coincide with the regions of large particle size and 
almost a constant N higher up in the flame. However in the lower part of 
the flame the maxima of ¢ and D,, coincide with the occurrence of a 
minimum of N. A _ tentative explanation for this seemingly paradoxical 


behaviour is given in the last chapter of the current work. 


A shown by Fig. 2.2b the particles start at a lower size down low in 
the flame and grow as they move upwards. The two main mechanisms of 
growth (agglomeration and surface growth) account for the steady increase of 


the particle size. The increase of $ with height may be due to particle 
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inception or surface growth mechanisms. The optical tests do not provide 
information about neither the existence nor the relative contribution of these 
processes to the soot loading. The béciat of smaller particles and large 
number concentration moves closer to the center line with increasing height. 
The particles grow during this stage and their number concentration decreases 
probably through agglomeration. Particle oxidation processes finally cause the 
cessation of the growth mechanisms as the particles move toward the upper 
portions of the flame. As we see the particle sizes in the intermediate and 
higher portions of the flame grow larger than 100 nm. This implies that the 
particles do not belong to the Rayleigh regime thus indicating that the profiles 
there are not as reliable. However low in the flame the particle sizes are 
clearly in the Rayleigh limit. A test of the applicability of Rayleigh theory is 


performed by measuring R. 


ij? defined as the ratio of the scattering coefficients 


for two different angles 8, and 8; i.e. 


Qyy (9) 
Sa ee 
HosvensOlay OF) 


When D/ S10 it should be Rj 71. A series of values of Rij for *thrée “different 
scattering angles (8,=45°, 8,=90° and 8,=135°) at the position of maximum 
volume fraction is given in the following table reproduced from ref. 84. 


TABLE 21 


Experimental Light Scattering Ratio Measurements * 


* From ref. 84 
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It is clear that low in the flame R;;=1, while higher up R;,>1 thus invalidating 


the adoption of Rayleigh theory. 


The depolarization ratio p, is defined as the ratio of the energy flux 
scattered in the horizontal plane (Q,,) to that in the vertical plane (Q,,) when 
the incident light is vertically polarized. The depolarization ratio is zero if 
the scatterers do not exhibit any material and/or geometrical anisotropy. From 
the optical measurements we can construct the following table 

TABLE wen 


Experimental Light Scattering Coefficient and Depolarization Ratio 
Measurements 


1.8x1073 


As we see, at a height z=30 mm there is a drastic decrease of p, from r, to 


fore (LactOr Olen oy: The value of p, at r<r, is” comparable to (that 
corresponding to the annular region. This behaviour of p, is discussed in 
Section 3.6. 


The optical method also provided data on the fluoresence cross 
section Q, at a wavelength of 514.5 nm, using a laser line at 488.0 nm as the 


excitation source. The fluorescence profiles were found to be broader than 


-14- 


the Q., profiles, extending closer to the center line. The maximum number 
concentration in the central core coincides with the maximum fluorescence 
level. It was also found that low in the flame (z<S5mm) the Q,, signal was 
due almost entirely to fluorescence. The fluorescence signal shifts towards 
the center line with increasing height until the entire central part of the flame 
exhibits fluorescence. A very similar fluorescence has been observed in 


premixed flames of various fuels*® as well as in other diffusion flames®’, and 


13,75 However 


has been attributed to various polycyclic aromatic hydrocarbons 
the existence of other fluorescing molecules (polyenes and polyynes according 
to Smyth et al.%°) which seem to be unconnected with the soot formation 


processes can not be excluded. 


A study of the velocity and temperature fields of the ethene flame 
was performed by Santoro et al®® The axial (v,) and radial (v,) components 
of velocity are shown in Figure 2.3. For heights below z=20 mm both v. and 
v, exhibit a maximum in the annular region of the flame. The strong 
entrainment of air inside the flame at these heights is indicated by the 
substantial radial velocity component. Mitchell et al.’? reported similar 
behavior in the lower portions of the axisymmetric laminar methane diffusion 


flame they studied. A strong acceleration of the flow along the z-axis is 


apparent throughout the flame showing the strong effect of buoyancy forces. 


Temperature profiles were obtained using fine wire thermocouples. The 
uncorrected temperatures are shown in Figure 2.4 for various heights. An 
estimate of the radiation corrections for the thermocouple measurements in a 
methane diffusion flame is given in ref. 90. This correction is of the order of 


5% and is not expected® to be very different for the ethylene flame studied. 
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Comparing the profiles of ¢ and T we see that the temperature maxima 
correspond to radial positions larger than r,. The radial distance between the 


positions of 9 and Tae varies from | to 2 mm. 


max x 


According to Mitchell et al.7? the luminous flame surface shown in 
Figure 2.5 marks the fuel rich side of the reaction zone, which is defined as 
the zone corresponding to the maximum temperature. This reaction zone has 
a blue hue the intensity of which decreases with increasing height above the 
burner mouth (not visible above z=5 mm). The fuel lean side of the reaction 
zone is dominated by the burnout of CO and represents a surface through 
which oxygen diffuses to the fuel rich side. The minimum concentrations of 
oxygen and the fuel exist on the reaction zone (T_.,), where the products of 
the combustion exhibit their maximum concentrations. 


Santoro et al.% 


used the particle velocity measurements to determine 
particle paths (trajectories) inside the flame. The™ particle trajectories@are 
influenced by thermophoresis because of the presence of high temperature 
gradients inside the flame (as high as 400 K/mm). These gradients give rise 
to thermophoretic velocities up to 4 cm/sec which tend to transport the 
particles into the interior of the flame. However when compared to the 
convective velocities, the thermophoretic velocities are much lower in 
magnitude higher in the flame and are more important at lower axial 
locations. The particle paths calculated for the ethene flame are shown in 
Figure 2.6A. Using the calculated path lines. and the velocity field 


1 


measurements, Santoro®! related each position in the flame to a residence time. 


The exit plane of the fuel nozzle was used as the initial particle location. 
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The residence time is plotted in Figure 2.6B as a function of the two 
coordinates used (z,r). This figure was used to transform spatial dependence of 
various variables such as 4, D,, and N, to residence time dependence. Figure 
2.7. shows soot volume gracue particle size D,, and particle number 
concentration N, corresponding to the region of maximum @ as functions of 


the residence time. Only the curves labeled 3.85 cc/s are the ones related to 


our current work. These profiles were calculated by Santoro using Mie 


4 


scattering analysis for a self-preserving distribution of polydisperse spheres?°, 


which is less limited compared to Rayleigh theory but involves’ more 
computational difficulties. It is apparent from Figure 2.7b that the particles 
grow traversing through the low and intermediate regions of the flame. At 
the same time the number concentration decreases through agglomeration, up to 
the point where oxidation mechanisms take over. From then on, ¢ and D,, 
decrease, while N is observed to decrease as the particles move into the 
oxidation zone. “This behavior of the number concentration is different from 
what Kent and Wagner? found for a laminar ethylene diffusion flame based 
on a Wolfhard-Parker burner. They concluded that the number density at the 


region of maximum @ remains invariant with height. 


The information provided by the optical data is used in the third 
chapter of the current work in combination with the results obtained by our 
thermophoretic method, in order to examine the aerosol dynamics of the soot 
aerosol in the laminar axisymmetric ethene diffusion flame under 


investigation. 


=] 


CHAPTER 3 


EXPERIMENTAL 


A novel experimental apparatus was designed to probe an 
axisymmetric laminar ethene-air diffusion flame at atmospheric pressure. The 


83,84,86 and a wealth of 


same flame was extensively studied by optical methods 
data was produced as a result of these studies. The thermophoretic method 
presented here complements the optical method providing useful information on 
the morphology of the soot particles throughout the flame. The combination 
of the two methods can contribute to a qualitative understanding of the 
complicated phenomena associated with the soot formation processes. In the 
first part of this chapter the theoretical background is given for a thorough 
understanding of the thermophoretic mechanism that drives the soot particles to 


the sampling surface. In the latter part the apparatus used is described and 


the experimental procedure of the sampling method is presented in detail. 


3.1 Thermophoretic probes 


Small particles inside a nonisothermal environment are driven from 
the high to the low temperature regions. This phenomenon of particle transport 
down the temperature gradient was given the name thermophoresis, which 
literally means being carried by heat. Many particle sampling methods are 
based on the thermophoretic mechanism. More recently the role of 
thermophoretic deposition?***:79 has provided a quantitative understanding of 
the drift of particles across the boundary layer to a surface that is at a lower 
temperature than the surrounding gases. This principle has been used 
extensively *9»74:76,80,88,95 + collect soot particles from a flame environment on 


cold surfaces for subsequent analysis by electron microscopy. 
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The thermophoretic deposition is driven by the _ presence of a 
temperature gradient in the vicinity of a cold wall inside the flow field of a 
particle-laden gas. This gradient is readily established by introducing briefly 
into the hot flame gases a probe surface which is initially at room 
temperature. The probe exposure time T, should be long enough to capture a 
significant sample but short enough to present a cold surface to the flame-born 
particles. This cold surface serves a second important role which is that it 
freezes heterogeneous reactions of the particles that are already captured. This 
chemical freezing action inhibits changes in the soot morphology after the 


particles have impacted upon the cold surface. 


The een’ surfaces employed in this work are the two devices that 
are used to hold specimens in electron microscopy as well as modified versions 
of these materials. These devices are the widely used circular perforated grid 
of 3.05 mm diameter and the bulk specimen carrier, a thin strip of metal with 
dimensions 11.0x3.5x0.1 mm, used to examine bulk specimens’ through 
scanning electron microscopy (SEM). The more familiar porous grid, Figure 
3.la, is advantageous for transmission electron microscopy (TEM) work, 
especially when it supports a 200 A thick film of elemental carbon or silicon 
monoxide (SiO). These films are deposited on the grids by low pressure 
evaporative deposition and are commercially available in a wide variety of 
dimensions of the supporting grid. The carbon films are extremely stable under 
the electron beam, possess a fine background structure giving high-resolution 
capabilities, and have a composition closely related to that of the soot 
particles. Carbon films however are readily oxidized in the flame environment 


and therefore can be exposed only for very short time intervals. Films of SiO 
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are generally more durable but tend to produce a lower background contrast. 
Two different kinds of grids, supplied by Ernest F. Fullam Ines Latham NY, 
were used in the tests. These were the 200-mesh copper grids, made up of 40 
um bars forming 85 um square holes and the 400-mesh copper grids made up 


of 33 um bars forming 30 um square holes. 


The image quality produced by photographing soot particles captured 
on carbon substrates is very high and good particle statistics are potentially 
obtainable from the electron micrographs. A disadvantage of the porous grid is 
its inability to afford good spatial correspondence between the flame position 
examined and the electron micrograph coordinates. Nevertheless position 
correspondence to about 10 mm between flame and micrograph coordinates is 


achievable with this device. 


The bulk specimen carrier shown in Earine 3.1b, is a thin strip of 
metal which can be used to provide better position correspondence by aligning 
its edge with the flame axis as illustrated in Figure 3.2. This figure shows a 
schematic diagram of the probe inserted into the diffusion flame, and the 
burner based coordinate system applicable for our experiments. The bulk 
specimen carrier is commercially available by Philips Electronic Instruments 
Inc. (part no. PW6156/00). The material of the carrier is an alloy of copper, 
nickel and zinc. The original carrier was used initially to sample the particles 
from the interior of the flame. However the smooth surface of the probe was 
oxidized during the exposure to the hostile flame environment thus limiting the 
capabilities of the SEM observations. In order to provide a chemically inert 


surface to capture the particles, the carrier was coated with a gold layer of 
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nominal thickness of 200 A. The distance from the end of the Bar cick when 
placed in the microscope provides a reference to link the flame radial position 
and the micrograph coordinates to within about 0.2 mm. Magnifications up to 
50,000 with adequate resolution are obtainable in the SEM mode. A 
disadvantage of the use of the bulk specimen carrier is that it allows only 
SEM observations which provide poor image quality (low contrast) of the 


low-density soot particles and does not yield quantitative information on the 


particle morphology. 


Modifications of these probes have been made in order to combine the 
unique advantages that each affords. The first modification consists of 
replacing a sector of the specimen carrier with a 200-mesh nickel screen (type 
I probe in Figure 3.lc). The grid material is coated with a carbon film by 
evaporative deposition after being bonded to the modified carrier using epoxy 
cement. This probe offers good position correspondence between flame and 
micrograph coordinates as well as the high image quality afforded by TEM. 
The relatively large screen surface makes this probe fragile and more difficult 
to fabricate and coat uniformly. For this reason the carbon layer was 
supported by a thin formvar (polyvinylformal) film of 200 A nominal 
thickness. The formvar films have a tendency to shrink slightly when first 
irradiated but when stabilized with an ultrathin carbon coating, are more 
resilient than carbon alone and very stable under the electron beam. When it is 
imperative that only pure carbon is present in the substrate material, the 
formvar can be removed by immersion in chloroform. The carbon fine on the 
type I probes were manufactured by Ted Pella Inc., Tustin CA, and the 
formvar film was not removed since the probes had to sustain high 


acceleration forces during the sampling procedure. Comparisons of particles 
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collected on these probes and on pure carbon coated grids showed no influence 
of the formvar film on soot morphology. 

A second modification consists of a carrier through which a slot has 
been milled (type II probe in Figure 3.1d). Conventional coated grids are 
bonded to the carrier with epoxy cement to provide a durable and more easily 
constructed probe. This probe has a larger thermal inertia than the type I 
probe, is imexpensively fabricated from the commercially available 
carbon-coated circular grids and the bulk specimen carrier, and is highly 


durable. 


A third modification of the carrier is shown in Figure 3.le (type III 
probe). A strip of the original carrier was removed along the probe thus 
providing a sharp edge as a stagnation flow collector. This probe affords high 
spatial positioning but gives information only on the primary particle size. The 
state of agglomeration is not revealed by the edge probe because of the high 


rate of deposition at a stagnation point. 


A combination of the previous probes is shown in Figure 3.1f (type 
IV probe). A carbon coated grid was bonded to the carrier at a certain 
distance from the tip of the probe. This probe combines two different 
sampling surfaces that collect particles corresponding to a certain axial location 


(z) above the mouth of the burner. 


The modified versions of the specimen carrier could not be 
accommodated by the appropriate specimen holder used to perform SEM 


observations on the original probe. In order to use TEM to examine the 
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modified probes some simple alterations of the holder were made. The soot 
samples collected on the probes have been examined on a Philips 420 Scanning 
Transmission electron microscope operating in either the SEM or the TEM 
mode. This Me rnecn provides magnifications as high as 820,000 and possesses 
a point-to-point resolving power approaching 0.3 nm in the TEM mode. This 
capability is adequate to resolve the soot particles except in the very lowest 
portion of the flame where the smallest discernible size is determined by the 
instrument resolving power. Soot particles.of 1-2 nm would be discernible, 


but smaller particles could be indistinguishable from the grainy background 


structure of the carbon substrate. 


.2 Boundary layer studies of the flow field around the probe 


A numerical investigation of the flow field around the probe was 
performed using .a computer program called FLUENT, developed by Creare 
Inc., Hanover NH, solving the Navier-Stokes equations without any 
approximations. The problem considered is the two dimensional steady flow 
around a flat plate of finite dimensions, maintained at a constant temperature 
different than that of the surrounding fluid. This is a first approximation to 
the actual 3-D unsteady flow problem around the specimen carrier inserted 
into the flame. The thermodynamic properties of the combustion gases were 
approximated by those of nitrogen which was considered an ideal gas. The 
heat capacity dependence on temperature was expressed as a _ polynomial 
function of T. Data were obtained from the JANAF tables®®, while viscosity 
and thermal conductivity data were taken from Svehla9!. The free stream gas 
velocity and temperature were U=0.7 m/s and T.=1!1700 K_ respectively 


(characteristic®® of the flame environment at an axial location z=l10 mm) while 
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the probe temperature was taken T,=380 K from _ probe temperature 
measurements at the same Z. At this location we calculated a Reynolds 
number Re, =9 based on the streamwise dimension of the plate probe, 1.e. L=3.5 
mm. The numerical solution allows the full thermal boundary layer to be 
calculated. Figures 3.3A and 3.3B show the temperature and velocity contours 
around the probe. We define the thermal boundary layer by the isotherm 
T,, such that T,)/To=0.9. It is clear from Figure 3.3A that the thermal 
boundary layer thickness 6, at the low Reynolds numbers characteristic of the 
problem, is rather larger than expected (typical value §,~2.3 mm halfway along 
the probe, i. x= L/2). The calculation of 6, permits the particle trajectories 
to be calculated as well. The trajectory shown in Figure 3.3 is that 


corresponding to a particle ultimately captured at x=L/2. 


A time scale of importance, the quench time to is the time required 
for a soot particle to traverse from the hot flame gases to the cold probe 
surface through the thermal boundary layer. Figure 3.3C shows a typical 
particle P traversing the distance between two neighboring isotherms T, and 


1, The thermophoretic velocity V, is perpendicular to the fluid velocity 


n+l1° 


Unnsi giving a total velocity vector V toward the probe surface. The time 


required for the particle to move from P to Pp is given by 


At Soa (341) 


where (PM) is the distance traversed by the particle along the x-direction 


. . Uy 
during the motion from P to P . 
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The thermophoretic velocity driving the particles of all sizes is given by 


vT 
V, = -Kv (3.2) 
iF 


where v is the kinematic viscosity of the gas mixture, and T is the local ° 
temperature. The thermophoretic velocity coefficient K depends on the 


Knudsen number, defined by 


d 

Kn = aaa (3:3) 
where \ is the mean free path of the gas and R is the radius of a particle. 
According to Talbot et al.°?, the values of K are 0.30, 0.51, 0.54 and 0.55 for 
Knudsen numbers of 0.1, 1.0, 10. and 100. respectively. Thus for particles in 
the free molecular regime (Kn > 1) K = 0.55. For particles in the transition 
regime (Kn = 1) the value of K is comparable to that of particles in the free 
molecular regime. A typical value of \ for the ethylene flame studied is 


21,22 the 


within the range of 400 nm to 600 nm. From our previous studies 
largest soot aggregates were found to be as long as 300 nm. These particles 
belong to the transition regime but according to the above results, have almost 
the same thermophoretic velocity with the smaller particles which belong to the 
free molecular regime. This guarantees the absence of preferential deposition 


of particles of different sizes since they undergo similar theromophoretic 


transport inside the flame environment. 


The temperature gradient in Eq. (3.2) can be approximated by 


VT aay 1S 
T * (PN) T, 
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where (PN) is the distance between the isotherms T, and T,,, along the 
y-direction, T,=(T tT 4/2: and Vi Cas atic kinematic viscosity of the gas 
mixture ata! - Equation (3.1) was used to calculate the values of At a tt 


after graphical evaluation of (PM) from the contours in Figures 3.3A and 3.3B. 


The quench time is then given by 


to ) | (3.4) 


Using the values corresponding to the characteristics of the soot field at z=10 
mm we estimate from Eq. (3.4) bimphd (2 be approximately 10ms. The total 
distance along the x-direction, traversed by a typical particle collected at 
x=L/2 at the same height during the penetration through the thermal boundary 
layer, was estimated to be 6z=2.3 mm. This means that the particle 
ultimately captured at z=10 mm entered the thermal boundary layer caused by 
the plate probe at a height of z=7.7 mm. Thus the quenching effect of the 
probe, being felt by particles entering the thermal boudary layer, occurs at a 
relatively large distance (8z) upstream the flow field of the combustion gases 
low in the flame. However higher in the flame the thermal boundary layer is 
expected to be considerably thinner because of the higher velocities. A lower 
value of 6 would increase the thermophoretic velocities which means that 62 
Should be significantly reduced. For the type I and IV probes the value of 6z 
for particles captured at a particular axial location is expected to be very low 
because of the small values of 6 around the grid material. For the stagnation 


point probe (type III) 6z is expected to be minimal throughout the flame. 


Another time scale of importance is the disturbance time T,, which is 


defined as the time required for the flow field around the probe after its 


-26- 


insertion to achieve velocity values equal to 90% of the values corresponding 
to steady flow conditions. Under the assumption of instantaneous insertion of 
the probe inside the flame environment we calculated T,=3.5 ms, corresponding 
to conditions characteristic of the annular region at a height z=10 mm. The 
actual sampling times employed throughout the current work were longer than 
30 ms. This suggests that the effect of the transient disturbance introduced 


by the motion of the probe to the flow field around it should by very low. 


3.3 _Thermophoretic studies on soot particle transport 


The theoretical description of the thermophoretic transport of particles 
is provided by the analysis of the particle drift across a thermal boundary 
layer formed over a solid surface in a particle-laden, hot gas stream?4. 
Assuming spherical particles of uniform diameter D, in the free molecular 


regime, this analysis provides the particle flux to the cold wall as 


J 6Kv fy ie l dT 3 5) 
w _ os . 
. D,” Te ‘fo dy J 


wheres J, is the particle number flux, K is the thermophoretic velocity 
coefficient, v, is the kinematic viscosity of the host gas mixture at the outer 


edge of the gas thermal boundary layer, f,, is the particle volume fraction at 


e 
the outer edge of the gas thermal boundary layer, D, is the particle diameter, 
Tes se are the wall and carrier gas temperatures, « is the exponent appearing 


in the dependence of gas mixture thermal conductivity on temperature 


(kx=0.84), and (dT/dy) faz is the temperature gradient in the 
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cross stream direction at the wall. We approximate 


(3.6) 


We can then estimate the required exposure time T,, here defined as the time 
required for a 10% surface coverage of the bulk specimen carrier by soot 


particles, to be given by 


T= (3.7) 


Using in Eqs (3.6) and (3.7) values corresponding to the characteristics of the 


soot field at z=10 mm, ie., K=0.55 for free molecular flow*%?, 


V,=2.81 cm?/s 
corresponding to air at T,, f, =1.8x10°° cm?/cm? and T,=1700 K from ref. 86, 
T,=380 K from probe temperature measurements, D,=13nm from ref. 22 and 
§,=0.23 cm from the previous chapter, we estimate T, to be of the order of 32 
ms. In practice exposure times of that order or slightly longer are employed. 
It should be mentioned that Eq. (3.6) is based on the assumption that the 


aggregated particles are small compared with the mean free path of the gas 


molecules. This condition is fulfilled in most regions low in the flame. 
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3.4 Probe control system and burner 


The probe control system and the related equipment is shown in 
Figure 3.4. The precise spatial positioning of the probe within the flame is 
achieved by mounting the mechanical components on a rigid aluminum table. 
The burner (a) with the chimney (b) and the flow restrictor (c) are securely 
held by a post (d) with a vertical screw drive mechanism. The vertical 
adjustment of the probe with respect to the burner mouth is measured by a 
cathetometer, and changes in the height above the burner are measured by 
means of a dial indicator (e). The vertical position control of the 


thermophoretic probe is about 0.1 mm. 


The ideal probe control system would instantaneously insert the 
thermophoretic probe’ into the flame for a well defined, controllable time 
interval without causing any disturbance to the motion of the flame gate The 
actual probe control system has been designed and calibrated to move the 
probe with the maximum speed to a precisely defined position in the flame. 
The mechanical actuator for this system is a double-acting pneumatic cylinder 
(g) of 14-mm diameter and 25-mm stroke that is driven by an air pressure of 
4.5 atm through a pneumatic directional control valve (h). A_ translation 
potentiometer (j) supplied by SERVONIC, Costa Mesa CA, is mounted on the 
piston shaft and provides a precise record of the probe trajectory on the 
screen of a storage oscilloscope (k). The potentiometer is connected in series 
with a D.C. power supply (1.5V) and the oscilloscope. The motion of the shaft 
during the sampling procedure causes a variation of the resistance of the 
potentiometer, thus changing the voltage signal forwarded to the oscilloscope. 


The linear character of the circuit allows an excellent correspondence between 
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the position of the piston shaft and the time variation of the recorded 
oscilloscope signal. The reproducibility of the probe trajectory with a 
specified setting of the valve-timing circuit (i) proved to be very good. With 
this configuration we have been able to achieve transit times T, for the entry 
of the probe into the flame equal to 12 ms. The exposure time of the probe to 
the flame environment is controllable and can be reduced to as short as 30 ms. 
The entire probe trajectory is recorded by the linear transducer during each 
test. Figure 3.5 shows the probe trajectory and documents the brevity of the 
transit time as well as the precise spatial definition of the residence time. 
The burner is identical in design with that used previously®>*® 
for the extensive laser scattering and transmission measurements on the ethene 
coannular diffusion flame. A cross sectional drawing of the burner is shown in 
Figure 3.6. The burner brass tube (a) has a 11.1 mm id. and ‘is surrounded 
by the outer brass tube (b) of 101.6 mm id. An evenly-spaced, 14-mesh ceramic 
honeycomb (c) fills the annular space. The fuel tube extends 4 mm above the 
top of the burner. An airflow through the ceramic provides the coannular air 
stream which overventilates the diffusion flame. A mixing chamber (d) in the 
middle of the burner is filled with 3-mm_ glass beads (e), to eliminate 
nonuniformities in the air flow field. Just below the glass bead chamber there 
is a preliminary air chamber (f), supplying the air uniformly toward the upper 
parts of the burner arrangement. The two chambers are separated by the 
distributor plate (g) which is a blockage plate with the purpose to begin the 
process of achieving a uniform air flow. The beads are prevented- from 
penetrating through the holes of the distributor plate to the preliminary air 
chamber using a fine screen (h). The burner is shielded from laboratory air 


currents with a 420-mm long plastic transparent chimney (i) supporting a flow 
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restrictor plate (j) in order to reduce any air recirculation in the system. The 
flow restrictor consisted of a 143-mm diameter porcelain plate with a 23-mm 
diameter center hole and concentric rings of 5-mm diameter holes. The 
chimney could slide along the outside of the burner so only one hole provided 


access to all the axial locations of the diffusion flame. 


CP-grade ethylene (99.5% pure) and compressed dry air (maximum 
moisture 3 ppm) were directed to the burner system. Both gases were supplied 
by LINDE specialty gases. Throughout the tests reported herein, the volume 
flow rates of ethene and air duplicated the values used at NBS namely 3.85 
cm*/s and 713 cm?/s_ respectively. These flow rates were measured with 
rotameters, which were calibrated for each fuel using a soap bubble technique. 
Using these flow rates upon ignition, a flame of definite shape and height (88 
mm) resulted®>. This flame displays the bright yeilow incandescence 


characteristic of soot but releases no soot to the surroundings. 


A series of high-speed cine photographic studies were conducted to 
determine the influence of the probe motion on the flame. These tests showed 
that the probe enters the flame very smoothly but a flow induced by the 
probe support rod distorts the flame severely after the probe motion has 
ceased. This unsteadiness of the flame position destroys the _ position 
correspondence between flame and micrograph which has been a goal of our 
effort. To prevent the bending of the flame we have added a flow deflector 
plate, shown schematically in Figure 3.2, through which the probe enters the 
flame. The deflector plate remains stationary in space thus deflecting any air 


current induced by the motion of the probe support. Cine photography of the 
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probe equipped with the flow deflector shows the flame to be perfectly steady 
throughout the time period during which the probe is in motion. The flow 
deflector plate is always positioned parallel to the luminous envelope of the 
flame in order to minimize interference with the combustion processes. The 
distance of the plate from the luminous flame envelope is 4-5 mm _ which 
proved to cause no disturbance to the appearance of the flame. 
Macroscopically the insertion of the probe during the sampling procedure, was 
accompanied by the appearance of a thin long dark streak down stream the 
flow of the combustion gases. This streak appears as a discontinuity of the 
luminous zone of the flame, which starts from the position of the probe and 
continues to the tip of the flame. This phenomenon is probably a result of 
the quenching effect of the cold probe to the flame reactions. However the 
consistency of the results shows that the effect of this disturbance is not felt 


upstream by the soot particles ultimately captured on the probe. 


The maximum temperature of the carrier probe during the sampling 
has been measured by means of temperature-indicating liquids supplied by 
OMEGA Engineering Inc., Stamford CT. At z=10 mm above the burner 
mouth, where gas temperature has reached a local maximum® as high as 2000 
K, the carrier probe reaches a temperature of 385+10 K when the residence 
time is 65 ms. In an equal time period the grid material would assume a 
higher temperature since its thermal inertia is substantially lower. These tests 
also showed that the probe is subject to significant radiant heating prior to 
insertion into the flame unless a protective flow of cold air in the probe 
support annulus is sustained during this period. This cooling air flow is turned 
off about 1 sec before the extension of the probe. The high speed 


photographic studies proved this time to be adequately long for the flame to 
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recover from its influence. The transient burner mouth temperature was also 
measured immediately after ignition in order to determine the warm-up time 
of the system. A J-type thermocouple was used with a l-mm diameter bead 
which was adjusted on the outer surface of the fuel tube 1 mm below the 
burner lip. The readings were recorded every one minute after ignition (t=0) 
until steady state. This procedure showed that the system reaches a steady 
state after 6 min. The transient temperature of the surface of the flow 
deflector plate was also recorded when the deflector was positioned in place at 
an axial location z=20 mm above the burner mouth. These measurements 
verified the Breuione warm-up time of 6 min which was adopted throughout 


our studies. 


3.5 Experimental results and discussion 


The thermophoretic sampling technique Heel iped previously in detail, 
was used to study the morphology of soot aggregates collected from the 
coannular ethene diffusion flame. These aggregates consist of individual 
particles, called primary particles, of approximately spherical shape and rather 
uniform size. Our observations of the diffusion flame generated soot are 
consistent with the long established fact of the uniformity of the primary 
particle size per aggregate in the carbon black industry. According to a 


18,62,66 the distribution of the primary particles sizes is 


series of studies 
generally much narrower within a single aggregate than in the sample as a 
whole. The term aggregate is used from now on to depict a discrete, rigid 
entity composed of extensively coalesced primary particles fused together. The 


aggregates are highly rigid and stable under the high-energy electron beam 


(accelerating voltage of 80 KV) during the EM observations. 


=33= 


There is a difficulty in sampling the soot aerosol in the fuel-rich 
region of the flame. This region is surrounded by an annulus (~2mm thick) 
of high soot volume fraction® that peaks at r=r, (annular region). The value 
of r, decreases with increasing height, as shown in Figure 3.2. The probe 
surface sampling the regions corresponding to r<r, can be contaminated by 
particles from the region of high soot volume fraction during insertion. This 
suggests that observations of the morphology of particles collected at r<r, must 
be very carefully examined for the presence of soot originating from the 
annular region. . The soot material samples from rar, are free of such 


contamination. 


A vertical survey of primary particle size has been performed using 
various probes positioned inside the flame for exposure time T, varying from 
30 to 120 ms. . The longer probe exposure times were used at low regions of 
the flame, since the soot volume fraction there is considerably lower than that 
higher in the flame. However these long exposure times (~100 ms) were 
associated with partial oxidation of the carbon substrate. The micrographs 
examined were taken from the undamaged regions of the film. Analyzing the 
soot samples collected this way, we were able to verify the thermophoretic 
effect on particle deposition inside a flame environment. This effect is clearly 
shown in Figure 3.7, depicting soot aggregates collected on a carbon-coated grid 
from a height of z=30 mm. It is obvious from Figure 3.7 that the particle 
concentration is higher in the vicinity of the grid bar which is playing the 


role of the cold target attracting the particles. 


A series of TEM micrographs reproduced in Figure 3.8 shows typical 


particles which are present near r=r_, for the indicated values of height z. It 
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must be pointed out here that particle deposition occurs on both sides of the 
film. However, as Figure 3.8 shows, the low coverage of the film surface by 
soot aggregates, as well as the narrow depth of field at these microscopic 
scales, would clearly discern the majority of the particles deposited on opposite 
sides of the substrate. Based on an examination of these and similar 
micrographs we come to the following conclusions. Large aggregates dominate 
the particle morphology at the annular region of the flame at all heights. 
These aggregates consist of primary particles of rather uniform size at low and 
intermediate heights. The primary particle size strongly increases as the 
particles move toward higher regions of the flame (Figure 3.8A-D). These 
Observations clearly indicate the dominant effect of surface growth on soot 
morphology at low and intermediate heights of the flame. The initial effects 
of soot particles oxidation are apparent from a height of z=40 mm (Figure 3.8 
E) and ape: Large aggregates of considerably smaller but uniform primary 
particle size are apparent. As the particles move toward the highest regions 
of the flame, these large clusters dominate showing a smaller primary particle 
size and a higher degree of agglomeration (Figure 3.8F). This dramatic 
decrease of the primary particle size with increasing height above the burner 
shows the effect of the oxidative mechanisms on the particles as they move 


toward the tip of flame. 


A detailed examination of the micrographs of Figure 3.8 provides 
some information on the contributions of agglomeration and surface growth on 
the soot growth mechanisms in our diffusion flame system. Low in the flame 
(z<10 mm) there exist intense particle inception and intense agglomeration. This 
causes a rapid broadening of the size distribution and a reduction of the 


number concentration of aggregates immediately after the cessation of the 
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particle inception pulse. At intermediate heights of the flame the number 
concentration has decreased through agglomeration (Figure 3.9A), and surface 
growth (Figure 3.9B) becomes the dominant phenomenon. Heterogeneous 
reactions result in material deposition on the exposed surfaces of the primary 
particles that make up the aggregates, thus retaining their nearly spherical 
form while growing (Figure 3.9C). The uniformity of the primary particle 
size is thus maintained, since all the particles forming an aggregate are 
exposed to an environment that is locally uniform within the scale of the size 
of the aggregate. This growth mechanism is compatible with the series of 
micrographs of Figure 3.8 and it also agrees with Tesner’s observations on 


thermal decomposition of hydrocarbons on carbon black aggregates®. 


The above observations were repeated with identical results, by 
performing the same: measurements with the use of SiO-coated grids. There 
was no obvious dependence of primary particle size on the substrate material. 
Measurements were also performed with various exposure times (30-90 ms) of 
the carbon coated grids within the flame at a specified height of z=20 mm. 
The grids were inserted vertically inside the flame in order to minimize the 
distrubance introduced by the probe to the flow field of the combustion gases. 
Figures 3.10 and 3.11 show soot particles collected from the annular region of 
the diffusion flame at z=20 mm using various exposure times of the grids. No 
change of either the primary particle size or the degree of agglomeration with 
time is apparent. Some signs of oxidation of the carbon substrate are 
apparent in Figure 3.10C and D as well as in -Figure 3.11C and D. The 
oxidation is probably associated with the long exposure times (70 and 90 ms) 
which result in the appearance of nonuniformities on the substrate background. 


No signs of oxidation are observed for shorter probe exposure times (T,=30 and 
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50 ms). The above observations lead to the conclusion that changes of the 
primary particle size or the degree of agglomeration do not take place on the 
probe, provided that the probe surface coverage by the particles is of small 
order, 10% or less. When the coverage is higher no positive conclusion about 
the degree of agglomeration can be drawn since some of the aggregates may 


overlap. 


Further support for the conclusion of no agglomeration on the probe 
is provided by sampling a methane diffusion flame where soot concentrations 
are dramatically lower®® (factor of 10-15 in soot volume fraction). The flame 
probed was based on the same coannular burner described previously in detail. 
A methane flow rate of 5.49 cm*/sec was used while the ventilating coannular 
air flow rate was 784.5 cm/sec. These flow rates produce a flame of 55 
mm height which displays the yellow incandescence characteristic of soot. 
Figure 3.12 shows the particle morphology as a function of the radial position 
r for the methane diffusion flame at two different heights (z=15 mm and z=25 
mm). Very small aggregates and single particles are present low in the flame 
as shown in Figure 3.12C and D. arees aggregates are formed higher up as 
shown in Figure 3.12A and B. At both heights the micrograph on the left 
hand side corresponds to the annular region of the flame. It is also apparent 
from these figures that the particles outside of the annular region are far less 
agglomerated and display a substantially smaller primary particle size (Figure 
3.12B and D). The above observations prove that agglomeration on the probe 


does not occur in the case of the methane-air diffusion flame. 


Examining a series of electron micrographs at different heights on a 


SMI Micro-comp data acquisition system, we were able to produce quantitative 
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data on the primary particle size distribution as a function of height z above 
the burner, for particles collected from the annular region of the flame using 
various probes. With the use of a digitizing tablet the primary particle 
diameter of about 100 particles at each height was measured. A basic 
statistical analysis at each height, performed by the same system, produced the 
mean value D, and standard deviation o of the primary particle diameter as a 
function of z. The values of D, vs. z are plotted in Figure 3.13. These 
values of D, compare favorably with the sizes observed with the SEM 
observations on the soot samples collected with the bulk specimen carrier. The 
maximum size. shown in Figure 3.13 occurs at z=40 mm which coincides with 
the axial position of the maximum local soot volume fraction found by optical 
tests (see curve labeled "NS annulus" of Figure 7a in ref. 85 ). We 


consider this level to mark the lower boundary of the oxidation region. 


The evolution of the primary particle size distribution with height for 
particles transported on the annular region can be characterized by the ratio 
o/D,. From our measurements we obtain that 0/D,, = 0.22, 0.18, 0.15 and 0.73 
corresponding to heights z=10, 20, 30 and 40 mm respectively. The above 
values show that the primary particle size distribution narrows considerably ° 
with height above the burner mouth. This behavior clearly illustrates the 
dominant effect of surface growth reactions on soot morphology at low and 


intermediate heights of the flame. 


As we see from Figure 3.13, there are no data points at heights below 
z=10 mm. This region of the flame was probed at a height z=5 mm. The 
soot sample collected at this height showed an aggregated stage, and a 


maximum primary particle size almost identical to that at z=10 mm. However 
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the radial positions at which D, achieved a maximum (thermophoretically 
determined r.) showed no reproducibility, thus indicating that the lower 
portions of the flame cannot be sampled consistently using the thermophoretic _ 


method. 


Figure 3.13 provides information on the specific surface growth rate 
f(t) of the soot particles corresponding to ind The specific growth rate is 
defined as the mass deposition rate per unit of particulate surface area, and it 
is generally a function of the axial location z. Tire -aepencence on 92 
PTHEIates Bto a dependence on the residence time t of the particles through a 
series of particle trajectory measurements reported by Santoro et al.2° for the 
identical flame. The’specific’ surface’ growth rate is related” with the 
material density p,, and the t-gradient of the primary particle diameter 


through 


Pp 
{a eee 3.8 
: ) (3.8) 


This equation was derived assuming a monodisperse primary particle size 
distribution function at each height z. This assumption was validated from 
an extensive evaluation of the electron micrographs which showed that the 
primary particle el is almost uniform for each aggregate corresponding to a 
specific location of the flame. If v. the axial component of the flow 


velocity inside the flame at the annular region then 


d 
— = Vv (3.9) 
dt 
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Then Eq. (3.8) takes the form 


f= (3.10) 


The axial velocity profiles as a function of z and r are taken from ref. 86. 
The values of D,(2) are displayed in Figure 3.13 for the _ particles 
corresponding to the annular region. In order to calculate dD, /dz we used a 
second order polynomial fit of the values of D,(2) which was subsequently 
- differentiated algebraically. Usingarcd. (3.10) we calculated the values of f as 
a function of z for the particles in the annular region. The z-dependence 
was converted to t-dependence through the particle residence time profiles 


calculated by Santoro®! for the identical flame. 
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A soot density of 1.86 gr/cm*? was assumed in Eq. (3.10). Thus we retrieve 


the values of f that are given in the following table 


TABLE 3.1 


Selected Data for the Ethene Diffusion Flame. 


Z te Vv. 3 res. time t : f 

(mm) (mm) (cm/s) (sec) (gr/cm?s) 

5 5.10 45 0.020 : 

10 4.50 67 0.028 94x 61.00" 

15 4.00 90 0.034 EOrsexes Use 

20 3.60 110 0.039 Oe US 

30 3.10 140 0.047 Cer 10." 
35 - eb 0.050 6 og boo a 

40 2.70 165 0.054 : 

50 2.40 - 0.059 : 


* From refs. 85, 86 
>’ From ref. 86 
° From ref. 81 


The calculations performed above demonstrate the potential of the 
thermophoretic method to give a quantitative description of the surface growth 
rate of soot inside a flame. No assumptions were adopted about the size 
distribution of the soot aggregates or the regime these particles belong to. The 
above data for the values of f(t) are used later on as an input to a theoretical 
model, in order to interpret the information derived by both the 
thermophoretic and optical methods about the soot field in the ethene diffusion 


flame. 
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Some comparisons are now made of the specific surface growth rates 
of soot calculated in our diffusion flame with those reported recently, 
concerning the soot growth rates in various flame environments. Our. studies 
resulted in a maximum value of f=10.5x107-° gr/cm’s corresponding to an axial 


location z=15 mm above the burner mouth (t=0.034 sec). Harris and Weiner*®*9 


2; for a series of rich premixed 


reported a maximum value of f=8x107> gr/cm 
ethylene-air, flat flames. This value was calculated by them using purely 
optical observations of the soot partice field, adopting Rayleigh scattering 
theory and assuming a self-preserving size distribution of smooth spheres. Some 
uncertainty related to the choice of the refractive index of soot was estimated 
to be rather unimportant. Miller et al.” reported a similar value for a 
methane-air diffusion flame on a slot burner. The relatively higher value 


calculated for our system is compatible with the fact that methane has a lower 


sooting tendency than ethene under similar conditions. 


Surface growth of soot particles in premixed flames has_ been 
attributed to various species. Prado and Lahaye’? suggested that 
polyacetylenes or polyaromatic compounds are the principal growth species. 


Homann and Wagner?’ 


concluded that only the largest polyacetylenes are most 
likely the growth compounds. More recently Harris and Weiner®®?9 presented 
evidence that acetylene is the principal growth species in premixed ethylene 
flames, and reported first order surface growth rate constants. Even though 
the surface growth mechanism has been thoroughly investigated in premixed 


2,4,6,7,41,43,98 there is too little known about the same mechanism in a 


flames 
diffusion flame environment. However it is known®®*:79.99 that acetylene is a 


major intermediate species in a diffusion flame. In a recent study of the 
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acetylene production rates in a methane/air diffusion flame Miller et al.”° 
reported a significant loss mechanism for C,H, in regions where soot particles 
are detected in high concentrations. These observations suggest®® that acetylene 
may be the soot growth species in a diffusion flame environment. Under 
this assumption the surface growth rate can be expressed as 


f = K(c,H,] (3.11) 


where K is the rate constant for the surface reaction that converts C,H, into 
soot, [C,H,] the mole fraction of acetylene and j the reaction order. 
Unfortunately our data for f start from an axial location of z=10 mm while 
the mass spectrometric measurements reported by Miller®® include acetylene 
concentration profiles up to z=9 mm. It is not possible therefore to estimate 
the reaction order from a Inf vs. In{C,H,] plot for our diffusion flame. If 
we further assume first order surface growth reaction kinetics, as supported by 


41,94 we obtain 


previous studies 


fesoK [C, Hy] (3.12) 


We should stress that the surface growth process in diffusion flames 
may be fundamentally different from that in premixed flames. Some recent 
studies’? however showed no proof of such a difference. Using the mass 
‘spectrometric measurements by Miller®® for our diffusion flame system, we 
have [C,H,]=0.022 corresponding to r, at a height z=9 mm. Using the 


corresponding value for f at z= 10 mm from Table 3.1 we can calculate from 


Eq. (3.12) the surface growth reaction constant K=4x1073 gr/cm*s atm 
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corresponding to the annular region at a height of z=9 mm above the burner 
mouth where T=I1700K. These values compare favorably with Harris and 


Weiner’s estimates for "young soot" in a premixed ethene flame.*? 


The radial dependence of the soot morphological characteristics can be 
investigated using the various modified versions of the bulk specimen carrier. 
Examination of soot samples collected at three heights (z=10 mm, z=l15 mm and 
z=30 mm) using three different probes showed that the thermophoretically 
determined values of r, (radial position of maximum soot concentration) were 
consistently lower than those detected by the optical technique. The 
descrepancy in the value of r, was higher down low (z=10 mm), while it was 
substantially reduced higher in the flame (z=30 mm). The radial positions, up 
to which particles were detected at a specified height, were also lower than 
expected from the optical measurements. These observations lead us to the 
conclusion that the three dimensional character of the flow field around the 
probe critically invalidates the position correspondence between the original 
radial position of ‘the particles and the position at which they are ultimately 
captured on the probe. However our observations indicate that the radial 
dependence of the morphology of the soot particles maintains its original 
character for heights z>10 mm even though the steep concentration gradients 
are obviously obliterated by the insertion of the thermophoretic probe. For 
axial locations z<10 mm the flow field of the combustion gases was probably 
highly disturbed by tthe insersion of the thermophoretic probe. This 
disturbance, combined with the sensitive chemical character of the flame at 
these heights, created a low’ reproducibility of the thermophoretically 
determined values of r, in the lowest portion of the flame. This discussion 


indicates that radial changes of the soot morphology detected by the 
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thermophoretic method show trends corresponding to the actual ones in the 
soot aerosol for heights z>l10 mm. We must emphasize though that the 
"particle profile shift" discussed above should not affect the size of the 


primary particles and its dependence on both radial and axial coordinates. 


At all three heights sampled with the modified probes we found that 
the primary particle size as well as the degree of agglomeration peak at r=r, 
(its corresponding lower value). This is demonstrated by Figure 3.14 showing 
particles collected from an axial location of z=30 mm and at three different 
radial positions. The probe used was the type I probe. An interesting 
observation can be made related to the particles collected from the fuel side 
of the annular region (r<r,). As Figure 3.14 shows there is a difference in 
the nature of the particles corresponding to r<r, from those corresponding to 
rr. The particles at r<r, possess a transparent or less dense appearance which 
“may be associated with a liquid phase. A tentative explanation that these 
particles are associated with the onset of oxidation of the carbon substrate is 
not reasonable, since the oxidation should start from larger radial positions 
where the temperatures are higher and the concentrations of oxidative species 
(oxygen etc.) considerably larger. It is clear from Figure 3.14B and C that 
oxidation of the carbon substrate at positions r2r, does not take place. The 
particles corresponding to r<r,were stable under the electron beam when 
captured on a carbon substrate. However when captured with the edge probe 
(type III) some motion was observed when the particles were irradiated by the 
beam. Figure 3.15 shows soot particles collected from three different radial 
positions using the edge probe at a height z=10 mm. The particles on the 


fuel side of r(r<r,) were further coated by a liquid phase deposit as displayed 
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by .PIgurels.) a, The above deposit can be associated with condensation 
processes of various species that exist in the flame environment. This 
condensation may be a result of the cooling of these species induced locally by 
the cold probe during the sampling procedure. The particles on r=r, (Figure 
3.15B) display a slightly larger primary size but no sign of the liquid deposit 
is apparent. Figure 3.16 shows particles collected from the fuel side of r, 
with the type IV probe from a height of z=30 mm. It is evident from Figure 
3.16A that the particles collected on the edge of the probe were further coated 
by a liquid phase deposit which is not as dense as in Figure 3.15A (z=10mm). 
The primary particle size. shown in Figure 3.16A compares favorably with 
that of the neighboring particles shown in Figure 3.16B, which were collected 


on the grid supported carbon substrate of the probe. 


As Figure 3.14A shows, the particles inside the annular region (r<r.) 
are of smaller primary size and display a smaller degree of agglomeration 
compared with the particles of the annular region (r=r,). However, some 
clusters of larger primary particle size were present, probably collected from 
the annular region as suggested by their morphology (compare Figures 3.14A 
and 3.14B). This fact does not allow for definite conclusions on the 


morphological characteristics of the particles at n<tpe 


Another important observation refers to the particles outside of the 
annular region (r>r,). For these particles the primary size and the state of 
agglomeration decreases drastically with increasing radius. This is demonstrated 
by Figure 3.17 showing particles collected from an axial location z=30 mm at 
three different radial positions (r2r.) using the type LV; probe. This 


observation was consistent with all probes and all sampled heights (see Figures 
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3.15 and 3.18 corresponding to z=10 and z=15 mm respectively) as well as with 
our previous SEM observations using the original bulk specimen carrier as the 


thermophoretic probe. 


Statistical analysis of the morphological data of soot Polen from 
the annular region performed by the SMI Micro-comp data acquisition system 
yielded the mean primary particle diameter D, at the sampled heights. These 
values are shown in Figure 3.13 and are in very good agreement with the 


mean sizes found by using the grids. 


A detailed comparison of the particle sizes found through the 
thermophoretic method with those of the optical method®® requires quantitative 
data on the aggregate morphology. Basic definitions of various terms used in 
carbon black image analysis are discussed: in ref. 1. Various methods*®*have 
been developed in order to analyze the aggregate size and shape of carbon 
black particles by electron microscopy. The agreement between these methods 
has not always been satisfactory. Medalia’s formulation on the morphology of 


61-63,65 was adopted in our present work in order to acquire 


aggregates 
quantitative morphological data on the soot samples collected using the 
thermophoretic technique. If D, is the primary particle diameter then the 


projected area and the volume of a primary particle are given by 


nD, 
3 
aD, 
y 9 we (3.13b) 
P 6 


-47- 


According to Medalia’s computer simulated random floc studies, an aggregate 
of projected area A consists of NS primary particles given by 


1:15 
N, = (a/a,} | (3.14) 


This formula was found to be accurate enough for 2<N £136. The 


equivalent volume of the aggregate is then given by 
Vo a a NM | (Sabo) 
Considering the equivalent sphere of diameter D, and volume V, we obtain 
D.=D.N (3.16) 


A series of micrographs of particles collected from a height of z=15 
mm was analyzed in order to obtain data on the size distribution of the soot 
aggregates. The type I probe was used to thermophoretically capture the 
particles. Two sets of micrographs corresponding to two different radial 
positions (r=r, and r>r.) were prepared for the image analysis. Each set 
consisted of about 10 photographic prints with approximately 90 numbered 
aggregates. Generally each particle that appeared continuous was treated as a 
single aggregate. However the low areal coverage of the prints reduced 
substantially the problems from overlapping. Aggregates that were cut off by 
the border of the prints were not measured. Each one of the particles was 
individually traced using the light pen of a Dapple Image Analysis system in 
order to compute its projected area A. The previously measured Be. was 
used to calculate A, from Eq. (3.13a) and then Ms from Eq. (3.14). Finally 
Eqs. (3.15) and (3.16) were used to calculate the volume of the aggregate and 


the equivalent sphere diameter. Two typical micrographs, one corresponding 
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to the annular region and the other to the outside region of r, are reproduced 
in Figure 3.19A and 3.19B respectively. The primary particle diameters were 
found to be ®h = 22.1 nm at r = 3.05 mm =r, and D, = 1.3 nm at r = 3.73 
mm. Figure 3.19C shows a typical micrograph of soot particles while . Figure 


3.19D is the equivalent digitized picture created by the Dapple system on a 


video image display. 


In addition to the parameters described so far some other quantities 
were also calculated; the perimeter p of each aggregate, a shape factor defined 
by f, = 4nA/p? showing the proximity of the projection of the aggregate to 
the circular geometry (f, = 1 for a circle), and the maximum length of the 
aggregate passing through the center of mass of the projected area of the 
particle. A test of the reliability of the tracing method was performed using 
two different operators to analyze a selection of a few aggregates. For 
intermediate and large size aggregates the measurements for the projected area 
A were well within 10% of each other. For the smallest particles (~5% of 
the total number of particles measured) the values of A differed up to 50%. 
These Hiererencct are attributed to the poor definition of the smallest particles 
but are expected to contribute insignificantly to the evaluation of the 
generalized moment ratios of the PSDF. Generally the linear quantitites 
measured by the two operators compared very favorably, while area 


measurements were within the experimental error requirements. 


The procedure described above provides the distribution of the 
equivalent sphere diameters D,. In the following discussion the notation D. is 
used to denote the values assumed by the equivalent sphere diameter D, of the 


soot aggregates. For a general discrete distribution with k size classes and 
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n, particles per class, the geometric mean diameter D, canses be 


calculated by 


InD, I (3.17) 


7. = (3.18) 


For the data obtained through the image analysis system n=l. So Eqs. (3.17) 


and (3.18) simplify to 


Peek aye (3.19a) 


(3.19b) 


where k now is the number of particles measured. 
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Using the same formulation, the generalized moment ratios D,, and D,, can be 


calculated using the following relations 


3 l k 3 
Di = T= D, (3.20) 
i= 
k 6 
yD) 
3 ; ; 
i=1 
eg ieee (3.20b) 
k 3 
ep 


These two equations provide the tool for comparisons between the soot 
morphological data obtained using the thermophoretic method and the data 
acquired through the optical studies. Using the quantitative image analysis 


results we can construct the following table 


TABLE 3.2 


Radial Dependence of Soot Morphological Data for the Ethene 
Diffusion Flame 


D,)(nm) 


(3.20a) 


D,3(nm) 


(3.20b) 


hike 


Santoro. et using Mie theory, calculated a value D,,2110 nm 


corresponding to particles of the annular region (r=4.0 mm = r,) at a height 
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z=15 mm. The use of Rayleight theory produced®® a value of D,,=78 nm. 
These values are in fair agreement with that calculated through the 
thermophoretic technique. As discussed previously the radial position of 
the annular region detected by the optical method is consistently higher than 
that determined by the thermophoretic technique. The calculation of D,, 
from the optical method required an assumption for the value of the complex 
refractive index of the soot particles. The assumed value of m=1.57 - 0.561 
proved to be very reasonable. The number concentration of the soot aerosol 
can be derived from the optical data using an approximation of ‘he value 
assumed by the- ratio fy=(Dg/D,,)*- Santoro et al.2® adopted the value 
corresponding to a self-preserving form of the distribution function which is 
equal to 2.00. As we see from our measurements fy has a value of 1.84 and 


1.93 at the two radial positions analyzed. 


Dobbins and Mulholland?* modeled the soot aerosol using the 
assumption that all the gas to particle conversion is due to particle inception, 
and particle growth occurs only by coagulation. Under these assumptions the 
replenishment of incipient particles coexisting with coagulative growth results 
in a wide distribution function and very high values of fy. This 
formulation, which is not inconsistent with the optical observations, is clearly 
shown to be incorrect by the recent morphological data obtained using the 
thermophoretic sampling technique. The dominant role of surface growth has 
the opposing effect or the variation of fy. Our observations as well as the 
numerical simulations (see Chapter 4) show that surface growth causes a 
reduction of the width of the distribution function and lower values of the 
ratio fy. In fact this quantity throughout the low portions of the flame is 


not greatly different from the asymptotic value that was chosen by Santoro et 
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al.85 and others in reducing the optical data. The applicability of this value, 
which is fortuitous because the assumptions for its validity are clearly not 
fulfilled, now appears to be a plausible assumption on which to reduce the 


4 


optical data. 


The mean values of the shape factor Ue calculated at the two 
positions analyzed, show that the particles are closer to a spherical shape away 
from the annular region. This observation is compatible with .our previous 
comment on the radial dependence of the degree of agglomeration of the soot 
particles at a specified axial location. Two more points of interest are; the 
constant value of De (number of primary particles per aggregate) with radial 


position and the relatively high values of the mean standard deviation 0,=0.4 


at both locations. 


Some comparisons of the actual equivalent sphere diameter 
distributions are now being made with the equivalent log-normal distributions, 
defined as those giving the same values of D,, and D,, with these calculated 
in Table 3.2. For these distributions D= 65nm, 0,= 0.26 at r.; D2 34nm, 
O= O27 eat ror © Figure 3.20 shows a comparison of the actual with the 
equivalent log-normal distributions at both radial positions analyzed. Generally 
the descrepancies are higher at the extreme values of the diameter range. 
However the differences at small diameters is not expected to affect the 


resulting values of Dz, and Dg. On the contrary slight differences at large 


diameters would result in a substantial change of the values of D,, and Dg. 


A computer experiment was performed in order to check the 


statistical reliability of the distributions calculated previously. Ninety 
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diameters were randomly generated in the equivalent sphere diameter range 


defined by aD Sande. These are the extreme values obtained by the 


n 
thermophoretic method for particles corresponding to r, at z=15 mm. Te 
ninety diameters were used to calculate the corresponding probability function 
values P(D). The equivalent log-normal probability function P was_ used, 
described by D,=65 nm and 0,=0.26. Nine size classes were defined in order 
to obtain the discrete distribution of the computer generated sample. The 
histogram of the distribution obtained is compared td the actual log-normal 
form (Figure 3.21A) corresponding to the distribution of a very large sample 
of diameters. We can clearly see that even though the diameters were 
distributed according to the log-normal form, the distribution obtained is very 
different. Another experiment was performed using a sample of 500 randomly 
generated diameters. In this case the distribution obtained is much closer 
to the log-normal form (Figure 3.21B), thus demonstrating the need of the 
morphological analysis of about 500 particles at each position in order to 
obtain statistically reliable values of the quantities measured. The above 
discussion illustrates the relatively low statistical reliability of the soot 
morphological data obtained so far. However even if the quantities 
calculated were not measured accurately, we expect the trends to be right and 
the comparisons with the optical data fairly reasonable. A more detailed 
comparison of the data obtained by both methods is needed, in order to 
evaluate the soot particle morphological characteristics in the ethene diffusion 
flame. However this comparison requires extensive work which is not 


pursued at the present time. 
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3.6 A hypothetical description of the physical aspects of soot formation and 
transport in the laminar axisymmetric ethene diffusion flame. 


The data acquired through the optical and thermophoretic methods 
provide extensive information on the soot field inside the laminar ethene 
diffusion flame under investigation. These observations are now combined 
and evaluated in order to interpret the mechanisms of soot formation and 


transport in the lower and intermediate heights of the flame. 


During the discussion about the appearance of the ethene flame, we 
mentioned that the luminous yellow regions are surrounded on the air side by 
a thin blue hue which locates the main reaction zone. This zone, which is 
visible low in the flame but becomes indistinct higher up, contains the surface 
of maximum temperature that is named "flame sheet". The reaction zone also 


72,90 with the location of maximum concentration of the combustion 


coincides 
- products (CO,, HO etc.). The fuel side of the flame sheet defines the fuel 
rich region, while the air side is known as the fuel lean region. Figure 3922 
shows the reaction zone surrounding the flame regions where soot particles are 


detected. The locus of maximum soot volume fraction $,., (annular region) 


is also shown, corresponding to radial positions on the fuel rich side of the 


flame sheet. Particle paths as calculated by Santoro et al.8® are displayed 
throughout the flame. In the upperivportion -ofs, the. ft lamée{(z>20r! mm). the 
line of corresponds to a particle path which probably is nearly coincident 


max 
with the streamlines of the flow field. It is apparent from Figure 3.22 that 
the particle trajectories penetrate the soot annulus (r=r,) at axial locations z 
lower than 15 mm. From that height the particles are convected upwards 
following trajectories parallel to r,. Thise behavior sof the. particle 


c 


trajectories is due to the strong entrainment of air and the high temperature 
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gradients in the lower portions of the flame, and has not been thoroughly 
investigated so far in connection with the particle formation processes. | 

Almost a quarter of a century ago, Wolfhard and Parker19°,101 
examined spectroscopically various diffusion flames on the slot burner, which 
now has their name. They found that in the ethylene and methane/air 
diffusion flames they examined, the sooting zone occurs some millimeters to 
the fuel rich side of the reaction zone. Haynes and Wagner‘**, and Santoro et 
al.85 in recent studies of laminar ethylene/air diffusion flames observed that 
the particle number concentration increases toward the reaction zone. The 
number concentration subsequently decreases toward the fuel zone while the 
particle size increases, and the soot volume fraction achieves a maximum on 


(Figure 3.22). However, Santoro et al.85 found that the number 


Lae | 
Ul 
— 


concentration goes through a minimum and then increases towards smaller 
‘radial positions (see N profile in Figure 3.22). This behavior of N cannot be 
attributed to the use of Rayleigh theory for the reduction of the optical data 
in ref. 85, since the particle sizes low in the flame are clearly in the Rayleigh 


limit. 


The high number concentration of particles close to the reaction zone 
occurs under a negligible soot loading (low 9). This can only happen when 
the particles are very small (d<2nm) which suggests that the particles are 
formed locally in large numbers on the fuel rich side of the reaction zone. 
The radical and ion concentrations are highest close to the position where the 
temperature maximum occurs. These species have been associated with the 
particle inception mechanisms and are considered®> to be _ potential soot 


precursors. Additionally at the lowest portions of the flame where the 
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particle trajectories intersect the soot annulus, there exists a maximum of the 


Q,, signal’? 


occuring at r=r.. If we accept that the particles are. formed 
outside tr. close to the reaction zone, then these particles are transported into 
the cooler, more fuel rich interior of the flame by the strong convective as 
well as thermophoretic forces. During this transport surface growth becomes 
important since the temperature and radical/ion concentrations become too low 
to create more particles. This period is also associated with intense 
coagulation because of the high particle concentrations. This mechanism may 
account for the minimum of the number concentration at the radial position 
where @$ displays a maximum. The maximum of the value of ¢ is caused by 
the surface growth mechanisms which become more intense as the particles 
move from the zone of their formation towards the interior of the flame. 


8 


These regions contain various species as acetylene® and other unsaturated as 


well as polyaromatic hydrocarbons”, which are considered to be the major 
growth compounds. Mass spectrometric measurements on an identical flame 


performed by Miller® 


, showed that the acetylene production rate peaks at r=r. 
implying that the surface growth mechanisms at these locations are very 
intense. The above locations are also characterized by lower temperatures 
and display visible fluorescence during the optical tests. According to Tesner??, 
surface growth can continue to increase the soot loading (¢) at temperatures 
much below those required for particle inception. The transport of the 
particles into regions of even lower temperatures, associated with a decline of 
the reactivity of the gases, causes a steady decline of the surface growth rate. 


The particles also lose their reactivity?®? by various mechanisms which are not 


well established at this point. 


The above discussion suggests that the first particles, called incipient 
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particles, are formed on the fuel rich side of the reaction zone in a region 
where the temperature is a few hundred degrees lower than that on the flame 


2 at the lower 


sheet. The Q,, profile-changes quoted by Santoro and Miller® 
portion of the flame indicate that the characteristic time for the particle 
inception processes is of the order of 1 ms. According to Harris et al.*?, who 
studied a premixed flat ethylene/argon/oxygen flame, the particle inception 
process lasts for about 2-3 ms in that flame. It is generally believed that 


this time scale characterizes the soot particle inception mechanism in a flame 


environment. 


The following discussion gives a tentative explanation of the variation 
of the primary particle size of the aggregates along the radial coordinate at a 
specific axial location z. This theory is consistent with the previously 
presented mechanism of soot particle inception in the lower portions of the 


ethylene/air diffusion flame under investigation. 


Particles on the fuel side of r, were found to have a smaller primary 
size than those of r., at all heights examined using the thermophoretic method 
(see Figure 3.14). These particles were created low in the flame (z<4 mm) 
where surface growth is not a significant factor. At these heights the Q., 
signal was very weak. According to Santoro and Miller®? the first axial 
location displaying the familiar Q\, annular structure was found to be at z=4 
mm. The particles created below this height were transported to the interior of 
the flame by the combustion gases (see Figure 3.22). These regions (r<r,) are 
also associated with lower surface growth rates as mentioned previously. The 
above mechanism explains the smaller primary particle size at r<r, compared to 


that at r=r, at all heights. This model can also account for the higher 
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number concentration on the fuel side of the annular region r<r, (see N 
profile in Figure 3.22). This point is discussed in chapter 4, where a bimodal 
integral solution of the aerosol dynamic equation is presented, and subsequently 


used to model the soot aerosol processes in the ethene diffusion flame. 


Particles present on the annular region at lower axial locations of the 
flame (4<z<10 mm) were formed outside r, at even lower heights. During 


their transport toward r., these particles grow through gas phase surface 


reactions. The particles on r=r, at z=10 mm travel in a richer® 


growth 
species environment and are relatively younger than those existing on the 
annular region at z=5 mm. The competition between the two factors could 
account for the observation of the thermophoretic method that the maximum 
primary particle size is almost constant in that height interval. From that point 
and on (z>15 mm), the particles are convected upwards along path lines that 
do not intersect the soot annulus. During this transport fee further grow 


through agglomeration and surface growth reactions until they enter the upper 


part of the flame (z=40 mm) where oxidation mechanisms dominate. 


Particles on the air side of r, have been formed outside r, lower 
in the flame. They were subsequently convected upward thus never having 
traveled through the region of maximum concentration of the growth species. 
The primary particle size is therefore smaller than that of the particles at r.. 
Another possible explanation for the smaller size as well as the lower degree 
of agglomeration of these particles is the presence of oxidative species in 
ample concentrations in a high temperature environment. ipnescursvecics 
(mainly O,) diffuse into the interior of the flame thus attacking the various 


particles and compounds in an environment which promotes the oxidation 
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mechanisms. 


The existence of smaller particles on the air side of r,, which was 
consistent at all heights using the thermophoretic technique (see Figures 3.14, 
3.15, 3.17 and 3.18), is also in agreement with the optical observations. AS 
shown previously in Table 2.2 for a height z=30 mm, the depolarization ratio 
p, drops by a factor of 3 from r, to r>r, According to the definition of p, 
this suggests that the scatterers (soot particles) exhibit less material and/or 


geometrical anisotropy toward the air side of r.. 


The above discussion attempts to combine the information provided by 
both methods for the ethene/air diffusion flame, in order to describe the 
physical aspects of soot formation and transport in that flame. However there 
are still steps that are not well established (i.e. mumber concentration anomaly 
at low heights), and more studies are required in order to gain a complete 
understanding of the physical processes associated with the soot formation 


mechanisms in diffusion flames. 
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CHAPTER 4 


THEORETICAL 


4.1 A bimodal integral solution of the aerosol dynamic equation 


The wide range of data related to the soot aerosol acquired through 
the optical and the thermophoretic methods was extensively evaluated in the 
previous chapter. In this chapter a mathematical model examining the aerosol 
dynamics is discussed when particle formation, coagulation and surface growth 
are present. The term coagulation is used to describe the process during 
which two spherical particles collide and form another spherical particle. This 
mechanism is used as a first approximation of the agglomeration processes 
inside the soot aerosol. The model accounts for a time varying bimodal 
representation of the aerosol particle size distribution function and provides a 
_theoretical framework for the interpretation of the experimental data. The 
bimodal solution is an approximate solution since it forces the size distribution 
function to be of a prescribed form (logarithmic normal). However the 
properties of this distribution are time dependent thus accounting for the 
evolution of the aerosol particle sizes. Upon completion of the model the 
effects of the various parameters on the development of the distribution 
function are investigated. Some conclusions about the soot aerosol are drawn 
contributing toward a better understanding of the mechanisms leading to the 
formation and growth of soot in diffusion flame systems. The model can be 
‘used to study the physical processes acting upon the particulate phase in 
various aerosol systems in which intense particle generation mechanisms are 


present. 
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The model is developed from the general aerosol dynamic equation by 
an integral formulation that includes the above physical/chemical mechanisms. 
Other models currently used to study the aerosol behavior are briefly discussed 
and extensively compared to this model. These comparisons show that the 
model provides a tool which could predict the effects of various parameters on 
the behavior of an _ aerosol, to assist the interpretation of experimental 
investigations. The development of the present model is discussed below in 


detail. 


4.1.1 Description of the model 


The time varying particle volume distribution function (PVDF) is the 
main quantity to characterize the growth of soot particles. A very 
common assumption in the past has been that the particles are of uniform size 
(monodisperse). However the existence of agglomeration strongly argues against 
this assumption and assures its invalidity. Another assumption has been the self 
preserving form*! of the PVDF in the case of an aerosol undergoing 
coagulation by Brownian motion in the free molecular’ regime. This 


self-preserving form of the _ distribution function (SPDF) was_ studied 


36 37 


extensively by Lai et al.°? Graham and Homer*®, and Graham and Robinson®’, 
These studies concluded that the SPDF is an asymptotic solution of the 
coagulation equation that is independent of the initial distribution of a time 
evolving aerosol. A necessary condition though for the validity of the 
self-preserving form is the constancy of the particle volume fraction @$ as 
discussed by Lai et al®”. This condition is violated when particle formation 


or surface growth are present as in the case of the soot aerosol in diffusion 


flames. 
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Numerous studies of samples of soot aggregates collected from 
various flames showed the particle volume distribution function to be 
approximately logarithmic normal or Gaussian (corresponding to a log-normal 


distribution with low standard deviation): 


Wersborg et al.% studying a flat acetylene-oxygen low pressure 
flame concluded that the PVDF is_ nearly Gaussian as_ long 
as most of the particles are spherical and it becomes log-normal 


when chain formation becomes predominant. 


In the pyrolytic decomposition of benzene diluted in nitrogen at 
temperatures around 1400K and atmospheric pressure, Lahaye et al.*° 
reported that for small particles (soot formation stage) the PVDF is 
Gaussian, while for large particles a log-normal distribution is 


observed. 


Chippett and Gray!? studying soot aggregates extracted from 
atmospheric pressure acetylene-air premixed flames have shown the 


sizes to be distributed according to the logarithmic-normal law. 


Bockhorn et al.° (1981) cCamé "to ine same —conciusion for the 
PVDF for particles collected from atmospheric pressure 
propane-oxygen premixed flames and _ propane-oxygen flames to 


which hydrogen or ammonia were added. 


Bochkorn et al.2 (1986) found the log-normal form to describe 


better the PVDF derived from EM studies of soot, collected from 
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premixed hydrocarbon-oxygen low pressure flat flames using various 


hydrocarbons (C3He, C>H>; C6 He). 


The above observations under various experimental environments led to the 


wide use of the log-normal form as the best approximation to the PVDF.°%-®.%° 


A continuous size distribution n(v,t) of an aerosol undergoing 
interparticle collisions in the presence of particle inception (formation) is 
governed by the following nonlinear integro-differential equation 


eo 


aie n(v,t) = seit K(v,u) n(u,t) du + 


a K(u,v-u) n(u,t) n(v-u,t) du + S(v,t) (4.1) 
where n(v,t)dv is the number of particles in the particle volume range v to 
v+dv per unit volume. The first integral in the above equation represents 
the loss of particles from that volume range because of collisions with 
particles of any volume. The second integral represents the gain of particles 
in the specified volume interval as a result of collisions between smaller 
particles. The source term S(v,t) accounts for the particle formation 
mechanisms and specifies the number of particles created in that size range 
per unit volume and time. The collision frequency factor K(v,u) represents 
the number of collisions of particles of volume v with particles of volume u 
per unit concentration of each size and its form depends on the size of the 


particles. 
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The particle volume probability function (PVPF) is then defined by 


n (v,t) 


P(v,t)dv = 
(v,t) = 


where N is the number of pparticles of all sizes per unit volume of the 
aerosol. It is obvious that | P(v,t) dv™= I. 
) 


The «th moment of n(v,t) is defined by 


M,(t)= | v¥ n(v,t) dv = N | v¥ P(v,t)dv (4.2) 


(2) 
°o 


Multiplying Eq. (4.1) by vk and integrating over all sizes 


f(u,v,k) K(v,u) n(v,t) n(u,t) dvdu + | S(v,t)vK dv (4.3) 


ro) 


a1 ok ee 0 
where f(u,v,k) = Ot keel 
2uv ;k = 2 


From Eq. (4.2) we derive that Mj=N, and M,= ¢ the particle volume fraction (total 


volume of particles per unit volume of system). Another variable of interest 


-65< 


is the dimensionless second moment defined as 


The quantities defined above have a physical interpretation related to the 
optical properties of small absorbing particles in the Rayleigh regime. The 
extinction coefficient is proportional to $ (see Eq. 2.4) while the volumetric 
scattering cross section is proportional to M, (see Eq. 2.1). AS shown from 
the definition of fy, this quantity is used to arrive at the value of N from 


the combined measurements of extinction and scattering. 


A general analytic solution for the first three moments with an 
arbitrary initial condition, does not exist due to complexities involved in the 
integro-differential form of the corresponding equations and the particle size 
dependency of the collision frequency K(v,u). Simplifications imposed on the 
form of K(v,u) can lead to analytic solutions. However in the case of 
coagulation due to Brownian motion in the free molecular regime (particles 
much smaller than the mean free path of the molecules of the host gas) the 
complicated expression for K(v,u) imposes the search for numerical solutions. 


In this case 


K(v,u) = Key, g(u,v) (4.4) 


where Kem = 


~aG— 


1 1 . 1 
l | 
and g(u,v) = c 's + Vv i: | es by, haa | 2 


u v 


ites Ke 1.38 x109716 erg mole¢! kK"! (Boltzmann’s constant) 
P,,: particle density 


T: Temperature of the aerosol 


This weakly singular form of the collision frequency is integrable, thus 
introducing no irregularities in the numerical integration of the 


integro-differential equation (4.3). 


~The above form of the collision frequency factor is derived in the 
kinetic theory of gases for collision among particles that behave as rigid 
elastic spheres. This assumption is reasonable for the soot aerosol at the 
early stage of soot formation where the. particles may coalesce upon colliding. 
However at later stages the particles retain their structure and chainlike 
clusters are formed. These highly structured particles called agglomerates or 
aggregates further collide under an "agglomerate coagulation coefficient" 
generally different from the form given in Eq. (4.4). Some recent numerical 


1.73 examine the form of this coefficient 


Simulation studies by Mulholland et a 
using a relatively small number of particles (~500). The law describing the 
Brownian motion of the particles is given by Langevin’s equation and some 


fractal theory concepts are used in their analysis. 


A strong particle formation term in the presence of coagulation would 
obviously result in a bimodal distribution that produces a double humped 


particle size distribution function? These rest: peak sas "tue particle 


ay 


inception (formation) mode and its height is determined by the strength of the 
source term. Once particles of the particle inception mode collide with one 
another they create new particles forming the growth mode whose properties 


are time dependent. The overall PVDF can be decomposed as 


n(v,t)= n,(v,t) + n,(v,t) = N,P,(v,t) + N,P,(v,t) (4.5) 


the two terms corresponding to the inception and growth mode respectively. 


Substituting Eq. (4.5) into (4.3) using (4.4) 


S(v,t)vKdv (4.6) 
O 
where the collision integrals are defined by 
(k) : 
Ztm. = f(u,v,k) g(u,v) P,(v) P,(u) dvdu (i = 1,2) (4.7a) 
) O 
(k) 
Ztm, = f(u,v,k) g(u,v) P,(v) P,(u) dvdu (4.7b) 
re) ) 


Equation (4.7a) gives the collision integrals resulting from _ interparticle 


collisions between particles of the same mode (i = 1 for inception mode, 
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i = 2 for growth mode) while (4.7b) results from collisions between particles of 


different modes. 


The source term is decomposed SOVst) eS (1) P tv) 


where Pj(v) the time independent, source particle volume probability function 
and S,(t) gives the number of particles of all sizes created per unit 
volume and time. Then Eq. (4.6) becomes 


dM, Kem 2 (k) . 3 (k) (k) 
: 1.8 
N, Engi + N, “tm, + 2N,N, “fm, |. S,(t) |. Bey ay. (4, ) 


We now prescribe the form of the time dependent particle volume probabiltiy 


function (PVPF) to be log-normal for both modes 


2 
l 1] l In V/V; (t) 
P: (v,t) = ee EXD - a cme ent yas ’ (i = 1,2) (4.9) 
cee oreo. (t) Me 
cs 1 


where Vg the geometric mean particle volume and o, the mean standard 
i i 


deviation. Both v, and o, are time dependent with dependencies we wish to 
1 1 


g 
determine. We further assume that the source PVPF P.(v) is log-normal, 


described by v, and o,, both constant. 
° 


fe) 


& 


For a log-normal distribution described by Vo, Og = o,/3 we can 


show that 


: : pts 4.10 
Me a exp 9, (4.10) 
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For a bimodal distribution 


M, (t) = vk [a (v0 ~ ng (vt) [dv = My (t) + M2? (t) 
O 


where the notation Mi is adopted to denote the kth moment of the PVDF 


corresponding to the i*® mode (i = 1, 2). 


Using Eq. (4.10) for a bimodal distribution consisted of two log-normal modes 


we obtain 


M. N : a : kK N < re 4 
,(t) = 1¥g, exp aaa Og + 2", CXD ee of - (4.11) 


Generally Ve. and Og. (i = 1,2) are time dependent. The last equation gives 
1 1 ‘ 


the first three moments as 


N = N, 4 N, (4.12a) 
9 2 9 2 
1 2 2 2 2 2 
M, = M; +M; = Digg exp ec gs + NaN ae exp cae (4.12c) 


The quantities N., 9, and Mi denote the first three moments of the PVDF 


corresponding to the it? mode (i = 1,2). The assumption on the log-normal 
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behavior of P.(v) simplifies Eq. (4.8) with the use of Eq. (4.10) 


dM, Ken 2 (k) 2 tk) (k) 
re = Sayeed IN, “fm, - N2eim, - 2 Zim. ]- 
Sia ee 4.13 
(Ove exp = au (4.13) 
For k = 0 we obtain the equation for the particle number concentration 


dN Ken 2 (0) 2 (0) (0) 
rm 5 IN, <fm, NS 2tm, + eee ]+ S,(t) (4.14) 


The first three terms of the right hand side of the above equation describe 


the following collision mechanisms respectively 


I) @+ @ O 
@ particles of mode 1 

II) O+ 0} 
O particles of mode 2 

III) os Oe 6 


As we see mechanism II does not affect N, while mechanism III does not 


affect N.. Mechanism I affects both N, and N, in opposite directions. N, 


is reduced by two units while N, increases by one unit during a single 
(0) 


2 
collision. A loss of one unit is accounted for by a term K, /2 N, Z¢m,- 


2 (0) 
Then a gain of one unit is represented by the term - K, /2 N, afm, * The 


source term in Eq. (4.14) affects only N,. This reasoning leads to the 


a 


following equations for N, and N, using Eqs. (4.12a) and (4.14) 


1 (0) (0) 4 
= Ka N, Soom - oot Boyds) (4.i5a) 
dt 
dN, Ken 2 (0) 2 (0) 
— [Ne Zem, = Ni Zem, | (4.15b) 


Using Eq. (4.13) for k = 1 we obtain the equation for the particle volume 


fraction 


d¢ 


fe) 


9 2 
= S,(t) eae exp = Og (4.16) 


Applying Eq. (4.10) with k =1 for the first mode and differentiating with 


respect to time 


do, dN, 9 2 dae 
my Vg, exp |—— og (4.17) 


Equation (4.17) was derived assuming that the size distribution parameters Vo 
1 


and Og of the first mode are constant with respect to time. This assumption 
1 


seems reasonable as long as the particles are generated by a log-normal type 


source with time independent parameters v, and o, . Obviously v, =v 
85 & fe) 81 So 
and Og = 9 . Combining Eqs. (4.12b), (4.15a), (4.16) and (4.17) we obtain 
1 fo) 
d¢, 9 2 (0) (0) 
= = “KenNiVg, exp aa hy Saget, TEIN ioe (4.18) 
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From Eq. (4.13) for k = 2 we obtain the equation for the second moment 


een 2 (2) 2 (2) (2) 2 2 
7 Ny Zem_ + No Zpm, + 2NiNy Zem_ [+ Sy(t) vg exp 189, (4.19) 
dt D 1 « b oO 20 


Again applying Eq. (4.10) with k=2 for the first mode and differentiating 


with respect to time under the assumption v, = v and 0, = 9 
81 es 81 8. 
1 
bls er veaa lige [18 8 . (4.20) 
= Vee (oj ‘ 
dt dt = =a 


2 2 9) 0) 
“KimNiYg, exp [18 2, | zea, aN. iii | (4.21) 


Equations (4.15a), (4.15b), (4.18) and (4.21) form a closed system with the 
unknowns N,, No, $4; M3 at a specified t. The collision integrals for the 


k 
first mode Ak can be calculated explicitly using the known values of v, and 
1 


81 
Og - The collision integrals for the second mode Zee) are complicated 
1 2 
functions of Vo and Og which can be expressed as functions of the 
2 2 


unknowns through the assumption of the log-normal behaviour of the growth 
mode (see Appendix A). The system of the above equations can be integrated 


numerically using any set of arbitrary initial conditions. 
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Initial attempts to integrate the system of equations for N,, N,, 9, 
and M3 proved to be unsuccessful because of the wide range of thempeiane 
variables involved. Using CGS units a typical number concentration value is 
of the order of 107° -101? (cm-*) while a typical value for M3 can be as small 
as) 10520 stouslO came Cie) The extremely wide range of the unknowns 
created some difficulties during the numerical integration. A dimensionless 
formulation proved to be the remedy to this problem. Introducing a new 
variable N., defined as the total number of particles created by the source 


mechanisms over all times per unit volume, we define 


e 
fe) 
i 
Z 
Pn 
ga 
oa 
3 
a) 
| 
Q 
ga 
~ 
Ee 


N, 
xo (4.224) 
NE 
N, 
xX, = 
2 N, (4.22b) 
N 9 2 
$, A aL ae og. | 
X, = of poem ce a AE ei 4.2 
3 d, 9 : ( 2G 
BS “gi SNe a5 Op 


=J4= 


The dimensionless time is given by 


t 
t! = 4, 2 
18 ene) 


where T, is defined below. 


So far we did not make any assumption on the form of the time 
dependent source function S,(t). A Gaussian pulse is a first approximation 


to the particle formation rate by natural sources 


b =it_ay2 
S(t) = So. exp ; oe | miSeen{i(t) (4.23) 


where t_ is the time of maximum formation rate, and T, is the e! half width 
of the pulse. From the definition of the total number of particles created 


over all times per unit volume 
N_ = | S,(t) dt (4.24) 


where i. is the time of observation which can be in the time interval where 


S,(t) is still active. Using Eq. (4.23) in Eq. (4.24), we find 
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SH ee ) t ey OTN les 
pistes. SL) aoe) (4.25) 
T td 


2 ae 
where erf(y)-= c dx 
Von 


Oo 


Dividing Eqs. (4.15), (4.18) and (4.21) by N, , 6, and M) respectively and 


using Eqs. (4.22) and (4.23) we have the dimensionless form of the system 


(0) 0) SERRE * ey. 

= Kg NGt X12 pte eal a (4.268) 
at. ° 
dX, KynN 7 2 (0) 2. (0) 
Fe ae agin SGA (4.26b) 
dX, (0) (0) 
me =- KeaNot ss] Selene. + XZ fm, (4.26c) 
dx, KenNat 2 (2) 2. (2) (2) 

~ xiZem + XoZem, + 2X, X,Zem 
dt’ 2v 2exp [180, ¥ : 2 

1 


(0) (0) 
- K tm NoTX1 [X1Zem, + 2am (4.26d) 


This system does not include the terms describing the surface growth 
mechanism. The development of the terms corresponding to the gas phase 


surface reactions is presented below in detail. 
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4.1.2 Surface growth term 


Coagulation, the first important mechanism of particle growth, 
modifies the size distribution of an aerosol but it does not change the particle 
volume fraction. The other important mechanism of growth, caused by 
surface gas phase reactions, affects the size distribution and causes an increase 
in the particle volume fraction. This mechanism is referred to as surface 
growth. The process involves heterogeneous surface reactions which occur 
under high concentrations of the condensing species. The, rater of,4 these 
reactions depends on the exchange of matter and heat between a particle and 


the gas phase. 


As discussed in a previous chapter the thermophoretic technique 
measurements lead to the determination of the specific growth rate f(t) of the 
soot particles, defined as the rate of gas phase mass deposition on the surface 
of the particles per unit surface area. The derivation of the form of f(t) 
from the EM observations did not require any assumptions about the growth 
species or the growth mechanism. In the following analysis the form of f(t) 


is considered to be a- known function of time. 


We consider a _ polydisperse aerosol growing by surface growth 
reactions. Following Friedlander’s analysis®?? we define I(v,t) to be the 
particle current or number of particles per unit time and volume of gas 
passing the point v in the v-space. If n(v,t) the PVDF it can be shown that 

dn dl 


pee iy Seer 4.27 
at Ov \ 
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We now consider the time derivative of M, due to growth and using Eq. (4.27) 


eo eo eo 


On vat dl . 
= Bene vK n(v,t) dv|= |v¥ ae duien ce ene ed (4.28) 
dt ot g av 


& re) ) 0 


From Eq. (4.28) for k = 0 we have 


= -[I iy Tas 20 (4.29) 


In Eq. (4.29) we used the fact that there is no loss or gain of particles by 


growth at the two ends of the distribution. From Eq. (4.28) with k = 1 


werven| Vans [iv] =: “ I dv = I dv (4.30) 
O O 


The term in brackets was set equal to zero since there is no loss or gain of 
material by growth from the two ends of the distribution. Neglecting 
diffusion in comparison with particle migration in v-space, the particle current 


can be approximated by 


I (v,t) oe (4.31) 
v,t) =n 
dt 
For molecular bombardment in the free molecular regime*? 
dv 2K 
Perr. Cor Vv (4.32) 


where Cor the surface growth coefficient depending on the molecular 


volume and mass of the growth species, the partial pressure of the gas phase 
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Substituting .Eqs. .(4.31),and 


growth species as well as the temperature. 


(4.32) in (4.30) we get 


si (4.33) 


The total surface area of soot particles available for growth per unit volume 


of the aerosol is given by 
1 / 
n| n D? dv = (367)'> My. 


ce) 


Sy = (4.34) 


If P, the density of the particles and f the time dependent specific growth 


rate 
d f 
= S. | (4.35) 
t : P, 
Using Eqs. (4.34) and (4.35), we find 
do ig 
——| = (367) Ma/3 (4.36) 
8 °p 


Comparing Eqs. (4.33) and (4.36) we obtain 


1 
/3 
Cy, = (36m) (4.37) 
P 
From Eq. (4.28) with k = 2 
dM —o 
ci = - | ae dv = - [v7] + 2| Ivdv (4.38) 
dt z Ne) 
8 0 
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The term in brackets is again zero. Using Eqs. (4.31) and (4.32) 


(4.39) 


where Cor is given by Eq. (4.37). 


Equations (4.29), (4.36) and (4.39) describe the effect of surface growth on the 
first three moments of the PVDF. As we expect there is no change in the 
particle number concentration of the aerosol since deposition of material on 
the surface: of fhe existing particles should not influence N. However, surface 
growth causes changes in $ and M, since the size of the particles is affected 


by the gas phase deposition. 


From Eqs. (4.12a) and (4.29) 


dN | dN 


2 
| et (4.40) 


dt ‘, g g 


Since there is no eeiaer Neve between the influence of the surface growth 
mechanism and the number concentrations of the two modes, under the 
restriction expressed by Eq. (4.40), a major assumption is made at this point. 
The effect of the surface growth mechanism on the first (particle inception) 
mode is neglected. This assumption is reasonable for aerosol systems where 
particle inception and surface growth reactions do not coexist. Under this 
assumption and using the log-normal form of the growth mode Eqs. (4.29), 


(4.36) and (4.39) take the form 


———$| = 0 = (4.41a) 
dt g dt g 
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/ 2 
si WOE NAA IY Pex B12 4.41b 
dt ‘g A | a.) sete 
2 
ma? aC IN “Is ye (4.41c) 
= r 2V exp |—- Oo . c 
dt g g 85 y Bo 
‘ 1 i 
where Cor = (367] sees 
Pp 
The quantities v and o, for a log-normal particle size distribution are 


Bi Bi 


functions of the first three moments, given by 


2 l N iM; 
o = —— In| —,— |; i = 1,2 (4.42a) 
Si 9 o. 
1 
2 
> 
ns fe= p12 (4.42b) 
j if a 
NEG? Miptet 


Combining these equations for i = 1,2 and using Eqs. (4.22a) -(4.22d) we obtain 


2 2 1 ' X X, ae 
cae sn os aa 9 ue 2 ( : a) 
Xs 
2 
X, 
re = uae 7 Y (4.43b) 
2 2 
Xo X, 
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Starting from Eqs. (4.41), using (4.22) and (4.43) we obtain the dimensionless 
form of the contribution of surface growth to the first three moment equations 


of the PVDF 


dx, dX, 
= QO t= (4.44a) 
i] 
dt g dt g 
2 8 =I 
dX, af ly ly 1g 
= Cor 5 Koi Tye heed (4.44b) 
ati g 3 2 
vee exp | —— Og 
aul 5 5 
dx, ce lg lg lg 
= 2 Cor ; X, Xz, i (4.44c) 
dt' g /s 11 21) 
Vg exp Spee: 
1p, 
where Cor = (367] 5 
f°) 


We can now formulate the complete system of equations describing an 
aerosol system undergoing particle formation, coagulation, and surface growth. 


The system takes the form 


dX, (0) (0) Smax Ts 

— = K,,N.7, Xi X.Zem cee Tb + ——__—. f, (t) (4.45a) 
dt' 1 b N, 

dx, 1 2 (0) 2 (0) 

a KON re Ix, Zor “sO ZF (4.45) 
dt! 2 2 5 . 
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2 8 ty 

dX, (0) (0) PEL as lox Js nage nee 

ns KN T. x, [%Zem, + X2Z¢m, + (367) iat 

dt! oe 5 2 
eset 

dx, Ken IN ee 2 (2) 2 (2) (2) 

oe Ix, Zem, + X2 Zfm, + 2%1X%2Zem |- 

dt! a 

iv eXD [186, } 


(0 (0) 1), £T,X,/2X 
“KigNoT X1[XsZem, 7 EES + [ 


The collision integrals zai ) Aged] Weandae> Oe areasiven by Eqs. (4.7) and 
are functions of o, , vg. As we mentioned before 0, = o, and vy = v, 

i i 1 ° 1 So 
while one ‘5 can be expressed as functions of the unknown dimensionless 
moments X,, through Eqs. (4.43). Details about the evaluation of the collision 
integrals Zem\" are given in Appendix A. vee OOve a Systemes t 
equations can be integrated in time numerically using a standard integration 
package with any set of arbitrary initial conditions. We used a variable step 
fourth-fifth order Runge-Rutta method developed by Fehlberg?®. A schematic 


of the various processes involved in the derivation of the system of Eqs. 


(4.45) is given in Figure 4.1. 
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In order to study the relative importance of the particle formation 
and the surface growth mechanisms we define R, as * the “ratio. Gia 
integrated particle volume fraction due to surface growth over the integrated 


volume fraction due to particle formation, when both quantities refer to the 


same time. So 


9(t) . 
gS aa (4.46) 
3 o( 
(|, 
But t 
do 
o() | = * [a (4.47a) 
g i Pag 
O 
and | t 
dod 
o(t) | ~ dt (4.47b) 
f dt ‘¢ 
0 
From the previous discussion we can write 
aces | ene" | (4.48a) 
= 46a 
dt g dt g 
and 
Oe | oe 4.48b 
dt@r-aaem@dte ole ee 
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Using the dimensionless form of the variables and Eq. (4.17) we obtain 


dX, 
di! 
9 Ot ues 
Re(t’)= a (4.49) 


From Eq. (4.45a) 


dx, SHES) Ue, | men 
at f N, 
Substituting Eq. (4.50) in (4.49) and using (4.24) we obtain 
i] 
dX, 
R(t") = ieeh| ate (4.51) 
dt! ai¢ 
) 


The values of the integrand in the above equation are given by Eq. (4.44b) as 
a function of the dimensionless moments X, The integration is performed 


numerically using the trapezoidal rule. 


Another important parameter used to check the validity of the 
assumption of no growth on the particle inception mode is defined below. We 
consider the effect of coagulation and surface growth on the particle inception 


(formation) mode in order to compare their relative importance in the 


=p5 


formulation of the model. The volumetric loss rate of material from the 


particle formation mode due to surface growth reactions can be found from 


Eq. (4.36) and is given by 


*/s : 
?. = “Co pNiv exp [25 (4.52) 


where Cor is given by Eq. (4.37). 


The volumetric loss rate of material from the particle formation mode due to 


coagulation is obtained from Eqs. (4.17) and (4.15a) 


(0) (0) ae 
o. = K,,.N,(N, Ztm, + N2Zem, ) Vg exp |—— Og (4.53) 


The _ relative importance of the two mechanisms (coagulation and surface 


growth) on the particle inception mode can be studied by means of the ratio 


ie (4.54) 


which is proportional to the specific growth rate f (included in the 
expression for Cer) that is determined experimentally. Small values of R,. 
support the validity of our assumption. Higher values indicate that the 
effect of surface growth on the first mode cannot be neglected. During our 
numerical simulations the surface growth rate was introduced at a later time 
when the particle inception mechanisms had weakened substantially. This 


guarantees very low values of the ratio R,. 
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4.1.3 Transition from the bimodal to the monomodal integral solution 


As we mentioned before a strong particle formation term in the 
presence of coagulation would result in a bimodal PVDF. The absence of 


the source term would obviously result in a monomodal configuration?® 


mainly 
due to coagulative growth that results in the transfer of particles from the 
first to the second mode. This rapid decay of the particle inception mode in 
the absence of particle formation results in a monomodal form of the PVDF. 
The only mode describing the aerosol system is then the growth mode and its 


evolution is described by the following set of equations as derived from the 


system describing the bimodal solution when N, = 0 


dN, 1 2 (0) 
=— K,N, Zem, (4.55a) 
dt 2 
do, 2) 2 
= 3 
= CgpNovg_ > exp (29, ] (4.55b) 
dt 
dM? 2 (2) es 2S ee 
oes KemNo Zpm, + 2CgrNovg, ° exp(— 9, ] (4.55c) 
Is f 
where Cor = [367] ar (4.37) 
Pp 


The wide range of the variables involved in the above system of equations 


leads once again to a dimensionless formulation. We define 


y= (4.56a) 
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Yes (4.56b) 
>, 
M3 

Y; = (4.56c) 
Se 


where N, is the total number of particles created over all times per unit 
volume (see bimodal formulation) and $,, M) are defined in Eqs. (4.22c) and 


(4.22d). The log-normal form of the PVDF leads to the following relations 


2 2 ] YY; 
og. = oat + —In |—,- (4.5 7a) 
Y, 
2 
Y, 
Fe Ve. ee (4.57b) 
/o /9 
Y, Ys 
Using Eqs. (4.56) and (4.57) in (4.55) we obtain 
dY, 1 2 (0) 
zi 2 Kim No Ts 4 Zem (4.58a) 
dt; 
dY 6 ; : im 
2 gr 's lg lg c! 
= ———————————— Y, 9Y, | Yg (4.58b) 
dt! 


-88- 


s 


dY, Kats 2 (2) 


2 2 
2v,  exp|18 
ve exp| og 
9 
ones @s (4.58c) 


The transition from the bimodal configuration, described by Eqs. (4.45), to the 
monomodal configuration described by Eqs. (4.58) happens when the number 
concentration of the first mode (N,) is negligible compared to the number 
concentration of the second mode (N,). The current wate: GIN Sah are 
used as the initial values of Y,, Y, and Y, for the numerical integration 


of Eqs. (4.58). 


An important parameter in the interpretation of the optical data, 


called the generalized mean diameter or moment ratio®*, is defined by 


foo) 


Bix DP dD 
P-q rAmsrqis ot 


Dig (4.59) 


fe D4 dD 
0 


1 
6v 3 
where D is the diameter of a particle of volume “[p = [ae| | and P(D) 


is the particle size probability function. The quantity P(D) dD represents the 
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fraction of all particles contained in the interval D to D + dD and has the 


property 


vo) 


[po dp = | (4,60) 
O 


The particle volume fraction can be written as 


Tl $ 
d= = ND, | (4.61) 
NM, 
From Eq. (4.59), using the definition of fy = : we can derive 
p 
3 6M, 6fy 
Da See 4.62 
63 mo mN ) ( ) 
From Eqs. (4.61) and (4.62) we obtain 
3 
Des 
fy = |--— (4.63) 
Ds 


; E D/Dg | 
¥20 
PD) = ee ee eee (4.64) 
27 Og D 


where D, is the geometric mean diameter and Og the geometric mean standard 
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deviation. Using Eq. (4.64) it can be shown that 
2 
p°+ q : 
Dog = Dz exP 222 og | (4.65) 


Using Eqs. (4.63) and (4.65) we find that for a logarithmic normal distribution, 


fy is related to oO, by 


2 
fy = exp [9g (4.66) 
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Symbol 


K(v,u) 


n(v,t) 


P(D) 


TABLE 4.1 


Nomenclature used in Bimodal Solution 


Meaning 


surface growth coefficient 


particle diameter 


geometric mean diameter 


generalized mean diameter 


specific surface growth rate 


particle current 


Boltzmann’s constant 


collision frequency factor 


free molecular coefficient 


th moment of PVDF 


particle volume distribution function 


particle number concentration 


number of particles created over all 
times per unit volume 


particle diameter probability function 
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CGS _ Units 


cm 


cm 


cm 


em7! 


Svmbol 


P(v,t) 


TABLE 4.1 (Continued) 


Meaning 


particle volume probability function 


growth to particle formation volume 
fraction ratio 


growth to coagulation volumetric 
loss rate ratio 


particle formation rate 
time dependent source function 
maximum value of S, 


total surface area of particles.per unit 
volume of aerosol 


time 

temperature 

time of maximum particle formation rate 
particle volume 

geometric mean particle volume 

first free molecular collision integral 


second free molecular collision integral 
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CGS _ Units 


cm 


Sec 


sec 


TABLE 4.1 (Continued) 


Symbol Meaning CGS Unis 
Greek 

Be particle material density grcm7 
Og geometric standard deviation in terms 


of particle diameter 


ao geometric standard deviation in terms 
of particle volume 


6) e-! half width of a Gaussian source pulse sec 


o particle volume fraction ‘ cm? /cm? 


Subscripts (superscripts for M,) 


0) refers to the particles created by the source 
l refers to the particle inception mode 

2 | refers to the growth mode 

b refers to the interaction of the two modes 


Dimensionless variables 


fy = NM,/¢? ; dimensionless second moment 


Cia=ot/Te : normalized time 


=94— 


TABLE 4.1 (Continued) 


normalized quantities when both modes 
are present 


normalized quantities when only the growth 
mode is present 
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42 A_monomodal integral solution of the aerosol dynamic equation 


In our studies of the aérosol dynamic equation Wwe \cain@ useertie 
monomodal integral solution as another approximate method of solution. This 
method represents the long time solution but it has an inherent weakness when 
a particle inception term is present; it falsely suppresses the particle formation 
mode which can be very significant when particle inception is strong. However 
this solution can accommodate all three physical mechanisms involved in our 


aerosol dynamic studies: Particle inception, coagulation and surface growth. 


The equations describing the evolution of the first three moments of 
the log-normal PVDF can be derived following the analysis of the previous 
chapter. The particle formation and surface growth terms are identical with 
those used in the bimodal solution while the coagulation terms are those 
describing the growth mode of the bimodal solution when the inception mode 


vanishes. So 


dN 1 2 (0) 
ee mee | Zem + S,(t) (4.67a) 
dt 


dd. De ae */s 2 
— = S(t) Ve. exp |—— oe + Cor N vg exp (20, ) (4.67b) 


dt - 


=96= 


dM 1 2 (2) 2 2 
=>KiN Ztm + Syt)vg exp (180, ] + 


dt - 


2 


2 Geean "ls a 4.67c) 
er Nv exp 5 Og (4.67¢ 


N, 9, M, : the first three moments of the PVDF 


Og 5 ee geometric mean volume and mean standard deviation 
e : of the particles created by the source term 
Og, Vg : geometric mean volume and mean standard deviation 
of the PVDF 
S,(t) : Gaussian source term given by Eq. (4.23) 
Cor Surface growth coefficient given by Eq. (4.37) 


We use a dimensionless formulation of the above system defining the following 


variables 

N 
Y, = (4.68a) 

N 

ie} 

9 2 
’ N Vg €xp [— Og ] 
y.*. a . (4.68b) 
2 d 2 
ie} 
ip ay. exp — Bi } 
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a 2 
Nv, exp (180, ] 


(4.68c) 
2 


My eM exp e 8 Ps } 


where N, is the total number of particles created over all times per unit 


volume defined as 


N, = S,(t) dt 
) 
t 
Using the dimensionless time t’ = =r the system (4.67) takes the form 
8 
dY, 1 2 = (0) 
= —K.N.T.Y, Zpm + S(t) (4.69a) 
dt! “ 
2 8 : 
dY, Cor Te ly ly dp 
= S(t) + —_—_—_——__ Y,_ Y, Y;, (4.69b) 
Ve eexD to 
Bo 85 
dy. K,»N.T, 2 (2) 
= VA Tae Ge SES) 
dt! 2 2 
ee exp [1805 
-1 5 5 
2 Cor ie Ig lg lg 
re Yq aX (4.69c) 


where S(t) = S,(t) Tae 
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In the derivation of the above equations we used the following expressions of 


Og, Vg because of the log-normal form of the PVDF 
2 2 i es nee 
fe} = 09 iy 4.70a 
& Bo 9 2 \ 
Y, 
2 
Y. 
i Eg 9 ee 4.70b 
2 2 
Y, Ys 


The collision integrals rae (k=0,2) according to Eq. (4.7) are functions of Oy, 
Ve which are related to Y, through Eqs. (4.70). These integrals are calculated 
by various interpolation formulae presented in detail in Appendix A. The 
three coupled ordinary differential equations (4.69) are _ integrated 
numerically in time using a standard variable step fourth-fifth order 


Runge-Kutta technique, developed by Fehlberg?®. 
As in the case of the bimodal solution, we define R, asethne sratio sol 
the integrated particle volume fraction due to surface growth over the 


integrated volume fraction due to particle inception at a specified time t. 


Using dimensionless variables it can be shown that 


R(t')= | ——|} dt! (4.71) 


The values of the integrand in the above equation are given by the last term 


of the RHS of Eq. (4.69b). 
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TABLE 42 


Nomenclature used in Monomodal Solution 


Symbol Meaning ; CGS Units 
Cor surface growth coefficient cm s° 
"sf 
Ke free molecular coefficient cin 7 s7# 
M, second moment of the PVDF cm? 
N particle number concentration cm7? 
N, number of particles created over all cm73 


times per unit volume 


R surface growth to particle formation 
volume fraction ratio 


S, (t) time dependent source function cm7 35-1 
t time sec 
Ge geometric mean particle volume cm? 
(0) (2) i ‘yan 
Zem A Zem first and second free molecular collision cm. 4,5¢m 
integrals 
Greek 
Og geometric standard deviation in terms of 


particle diameter 


siuu= 


TABLE 4.2 (Continued) 


Symbol Meaning CG tintts 
7; e~! half width of a Gaussian source sec 
pulse 
) particle volume fraction cm? /cm? 


Subscripts (superscripts for M,) 


fe) eters to thes particles creatreaubyutne source 


Dimensionless variables 


tl=. = : normalized time 
8 
“ N 
, N, 
d 
= $ normalized moments of the PVDF 
ce) 
M, 
Y, = 
M2 


a Ob Ao 


4.3 Other models for simulating aerosol dynamics 


Another existing model studying the aerosol dynamic behavior 1s 
presented in this section in order to make detailed comparisons with the 
integral solutions analyzed in the previous chapters. The model, called 
MAEROS, was developed by Gelbard and Seinfeld**** for simulating the 
evolution of aerosol size and chemical composition distributions resulting from 
coagulation, intraparticle chemical reaction, gas to particle _conversion and - 
particle generation by a source mechanism. The technique is based upon 
dividing the entire particle size domain into m arbitrarily located sections. The 
particle masses of section j range from Vj.1 to Vv; where j = 1, m. Imposing 
conservation of mass for each of the s components for all the processes 
described above and approximating’ a form of the distribution function as a 
garere within each section, results in a set of m by : differential equations 
for the mass concentrations Qik of component k in section! } (j = Uz). 
Keowee gee S).> cL nIsSeset of equations is integrated in time using the 
fourth-fifth order Runge-Kutta method developed by Fehlberg?®. The available 
code is limited in section refinement such that v,; 2 2 v,, and 5 ¢ m ¢ 20. 
In order to make comparisons with the integral solution we considered only 
one component (s = |) and assumed no intraparticle chemical reaction and no 
surface growth. The coagulation coefficient was assumed to be that 
corresponding to free molecular coagulative growth (see Eq. 4.4). The Gaussian 
form of the particle formation rate was used (Eq. 4.23) to describe the source 


term generating particles within the first section. The particle mass 


distribution function n(v,t) is defined such that n(v,t) dv is the number 


en 


concentration of particles in the mass range [v, v + dv] at time t. Since 


there is only one chemical component we have 


i fac 
j 
Q; (t) = Vv n(v,t) dv (4.72) 
Vie4 
where Q; is the mass concentration of aerosol particles in section j (j = 1,2,. 
. m). It is obvious that the volume fraction over all sizes is given by 
m 
l 
ot) = — ) Q; (t) (4.73) 
Pp 
j= 1 


- where P, is the density of the aerosol particles. 


This technique considers only conservation of mass (volume). So in case we 
want to derive some information about the other moments of the distribution 
function we have the relate Q,(t) to n(v,t). Following Gelbard and Seinfeld’s** 
analysis we choose the diameter D as the size variable of interest related to 


the mass v through 


After a series of assumptions it can be shown that the particle mass 
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distribution function in section j is given by 


n(v,t) = —__»-___ (4.75) 


Then the number concentration in section j is given by 


Vis 

j 
: Q; 1 ] 

N/ (t) = 1 CGS 9 als eee St af dee (4.76) 
Vv. Vv 
+) 0 
In 
Mjat VR 


The number concentration over all sizes is given by 


m 


N(t) = ) NI (t) (4.77) 
jel 


The sectional second moment, as defined through the PVDF corresponding to 


section j, can be obtained as 


V. 
j 
P v2 
M3 (t) = n(v,t) dv (4.78) 
Pp 
eet 
Using Eq. (4.75) we have 
Q; (y i v1) 
M3 (t) = ————__ (4.79) 
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The second moment of the PVDF is then 


m 
M,(t) = ) M3 (t) (4.80) 
j=l 
Equations (4.73), (4.77) and (4.80) provide the variables used to compare 
MAEROS to the integral solutions. We have to emphasize that the MAEROS 
solution does not assume a form for the distribution function. However the 
fact that it preserves only the volume, introduces a degree of approximation in 
the evaluation of the various moments of the PVDF except from the particle 


volume fraction. 


A test related to the above approximation was performed choosing the 
volume as the size variable of interest, instead of the diameter D. For this 


case the sectional approximation of the distribution function assumes the form 
vivir (4.81) 


which substantially differs from Eq. (4.75) that was derived using the 
diameter D as the size variable of interest. These two forms of n(v,t), i.e. 
Eqs. (4.75) and (4.81), were used in an example describing free molecular 
coagulation in the presence of a Gaussian source pulse with the following 


values of the primary variables 


T = 2000K, t,, = 0.01 sec, T,= 0.005 sec, N,= 1.735x 10*9 cm™® 


= 3 = = . . 
peoe 2 g/cm”, on 10 nm and Se 0.0 


The values of N and M, differ by as much as 5% from each other for the 


time interval 0 <¢ t ¢ 0.15 sec. This example clearly illustrates the 
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dependence of the results on the choice of the sectional form of the size 


distribution function. 


A more severe test of the assumptions introduced to calculate N and 
M, is performed for the case of a constant coagulation coefficient where an 
exact solution can be obtained analytically?*, A comparison was made for a 
constant source term creating monodisperse particles over a specified time 
interval (rectangular source pulse). The two different forms of the sectional 
distribution function given by Eas. (4.75) and (4.81) were used. It proved 
that the sectional form of aes distribution function when the diameter js 
chosen as the variable of interest (Eq. 4.75) gave results in good agreement 
with the exact solution for the number concentration (within 4%), while the 
values of the second moment were systematically in error within 12%. 
However when the volume was chosen as the variable of interest the agreement 
was less favorable. This test resulted in the choice of Eq. (4.75) over Ea: 


(4.81) for the sectional representation of the distribution function. 


The above examples clearly indicate the approximation introduced by 
the sectional form of the distribution function and address the fact that the 


technique should be treated as another approximate solution. 


The sectional multicomponent aerosol model has been significantly 
expanded recently by Warren and Seinfeld9’. The new model, called ESMAP, 
incorporates the mechanisms of homogeneous:nucleation and condensation that 
are coupled to the vapor phase and conserve particle number concentration. 
The model allows examination of the balance between new particle formation 


and particle growth as well as coagulation and deposition mechanisms. 
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A comparative review by Seigneur et al.8? of mathematical models of 
aerosol dynamics was recently published. Three different approaches were 
considered based on continuous, discrete (sectional) and parametrized 
(log-normal) representations of the aerosol size distribution. TrPnew rests: OF 
this study provide useful information for the selection of an aerosol model, 
depending on the accuracy requirements and computational constraints 


associated with a specific application. 


TABLE 4.3 


Nomenclature used in MAEROS Solution 


Symbol Meaning CGS _ Units 
D particle diameter cm 
M} sectional second moment of the cm? 


PVDF for’ particles’in section "} 


M, second moment of the PVDF cm? 
n(v,t) particle mass distribution function gr-tcm~? 
N) particle number concentration in © cm? 
section j 

N particle number concentration over cm? 
all sections 

Q; mass concentration of particles in grcm7? 


section j 
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Svmbol 


Qik 


Greek 


TABLE 4.3 (Continued) 


Meaning 


mass concentration ofcomponent k 
in section j 


time 


particle mass 


particle material density - 


particle volume fraction 
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CGS _ Units 


grcm7? 


sec 


gr 


grcm~? 


cm? /cm? 


4.4 Comparisons of the models 


The two forms of the integral solution (monomodal and bimodal) are 
now subject to a detailed comparison with the MAEROS solution in order to 
investigate the reliability of the data obtained by both methods. A number 


of special cases is considered and the results are presented below. 


The initial value problem (no source term) with free molecular 
collision frequency is examined first. For this example only the monomodal 
form of the integral solution is used since there is no particle generation term. 
In this case both solutions are approximate and no positive statement can be 
made for the accuracy of either one, since no analytical solution is available. 


In the example considered, the time evolution of an aerosol undergoing free 


molecular coagulation is_ studied. The initial number concentration is 
PNeeteuslOVeicmirs dndethe! particles srof density Pp = 2 gr/cm® are initially 
monodisperse with D = 10 nm. The temperature of the aerosol is T = 2000K 


(flame-like conditions). Twenty sections geometrically spaced in _ particle 
diameter from 8.96 to 911 nm were used for the MAEROS solution. ve 
particles formed by the source term had a diameter corresponding to the 
average diameter of the first section which was D = 10 nm. The values of 
N obtained by both methods were within 2% while the values of M, were 


within 17% for the time interval 0 <¢ t ¢ 0.23 sec. 


Another test for the accuracy of the two methods examined, is performed by 


the asymptotic value of the dimensionless second moment fy defined by 


NM, 


? 
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The MAEROS solution gives fy = 2.66 while the monomodal integral solution 
gives fy = 2.27, both at t = 0.23 sec. These values are in error Dy 2e7eaae 
9% respectively from the value fy = 2.079 recommended by Graham and 
Robinson?” for a normalized form of the self-preserving distribution function 
obtained by Lai et al5”) As we mentioned previously, the self-preserving 
distribution function is the asymptotic form of the PVDF for an aerosol 
undergoing free molecular coagulation by Brownian motion. The asymptotic 


value for the monomodal integral solution?? 


is fye = 2.297, a 10.5% discrepancy 
from Graham and Robinson’s value. However the MAEROS solution fails to 
give a reliable asymptotic value for f, because of a finite upper limit loss 
due to rapid coagulative growth of the particles at relatively high 
concentrations. The rapid increase of the particle size causes the early arrival 
of a large number of particles to the upper end of the diameter range, soon 
producing an undesirable loss of particulate material from the end of the 
distribution. The finite upper limit loss results in the inability of the code 


to comply with the volume preservation requirement thus distorting the 


predictions of the model. 


A more thorough comparison of the integral solutions with MAEROS 
is performed below for the case of free molecular coagulation with a Gaussian 
source pulse. The evolution of an aerosol is studied under "flame-like" 
temperature conditions; T = 2000 K. Particles with poa= 2 gr/cm?3, Dg = 5nm 
- and Co6 = 0.1 are generated by a Gaussian source pulse with t, = 0.03 sec, 
T,= 0.01 sec, shown in Figure 4.2A. This source term creates a steady increase 
of @ during the time interval in which it is still active. The particle 


volume fraction maintains a constant level after the time where the source 


becomes negligible as shown in Figure 4.2B. The total number of particles 
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created over all times is N, = 4 x 10/3 cm-* and the time interval studied 0 < 
thin 0:25 sec: Once again the unavailability of an exact solution does not 
permit any positive statement about the accuracy of the methods. ThewoP 
times for this particular case on a Micro-VAX were; 14 sec for the monomodal 
integral solution, 460 sec for the bimodal, and 214 sec for the MAEROS 
solution. A series of figures provides the comparison of the three methods. 
In Figure 4.3 the number concentration and second moment profiles for the 
three methods are compared. AS (we see the agreement. ofethe number 
concentration profiles is rather good, except during the time interval where the 
source term decays rapidly. The bimodal solution is closer to the MAEROS 
solution as expected, since the monomodal solution suppresses the particle 
formation mode thus distorting the actual form of the distribution. Finally the 
MAEROS solution gives consistently higher values for the second moment and 
the ratio fy as shown in Figure 4.4. It ‘1s ,alsow clears from, .the. previous 
figures that both forms of the integral solution tend to converge at very late 


times. 


The MAEROS solution is unable to give an asymptotic value of any 
of the variables, since the long simulation times are associated with a 
substantial loss of particulate material from the upper end of the size 
distribution. This fact is demonstrated by Figure 4.5 which shows the time 
Bertin of the mass concentration histograms of the aerosol. The bimodal 
configuration is clear at a time where the source term is strong (t = 0.02 sec) 
as shown in Figure 4.5A. The _ distribution assumes the monomodal 
configuration (Figure 4.5B) immediately after the particle formation term is 


turned off. The mass concentration histogram at t = 0.16 sec shows that the 
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largest particles have already reached the upper limit of the diameter range 
used in the calculation (Figure 4.5C), In fact the loss of material becomes 
important at a later time (t = 0.20 sec) where a substantial decay in the values 
of the volume fraction is observed. In order to avoid the finite upper limit 
loss a new definition of the diameter sections is required corresponding to 


larger sizes. 


Figure 4.6 presents the particle size distribution function n(v,t) 
computed with the two forms of the integral solution at a time t = 0.03 sec, 
when the source term achieves a maximum. The bimodal solution gives a 
very wide growth mode (Og. = 0.517) due to intense coagulation between the 
particles. The monomodal solution gives an even wider PVDF corresponding to 


o, = 0.603. The bimodal character of the PVDF is clearly indicated in the 


& 


figure, from which one can see that a log-normal or self preserving form 
would be obviously inappropriate to represent this distribution. Frenklach and 


Harris?? 


came to the same conclusion using a discrete model to examine a 
problem with simultaneous nucleation, coagulation and surface growth -(see Fig. 


6 in ref. 29). 


An important point is brought up by the form of the PVDF shown in 
Figure 4.6 for the bimodal solution. The particles formed by the source 


mechanism are characterized by a geometric mean volume v,. and a standard 


go 
deviation Og, which define the lower limit of the particle sizes involved in the 
calculation. Strong coagulation during the time when the source term is 


active creates a rapid broadening of the growth mode (Figure 4.9) which at 


g2 


0.517. This broad distribution in the numerical calculation creates an 


the time corresponding to Figure 4.6 is characterized by a value of a 
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artifact with no natural meaning i.e. particles of substantially smaller size 
than Vo are present ine relatively high concentrations as depicted clearly by 
Figure 4.6. However the effect of jthosé small particles in the, integral 
evaluation of N is surprisingly low as illustrated in Figure 4.3A which shows 
very good agreement of the profiles for all three solutions. This indicates 
the significance of the right hand side of the distribution function in Figure 
4.6 (corresponding to higher values of the volume) which almost coincides for 


both forms of the integral solution. 


hes s.importancesy-of utheespresences ,of ~-the.« smaller ‘-particles» in) the 
numerical evaluation of the higher moments is expected to be much lower, 


since n is multiplied by the volume v or v? 


for the calculation of $ and M, 
respectively. This above numerical artifact of the bimodal solution is present 
only during the time when the source term is active as shown in Figure 4.7 
which shows the PVPF of the growth mode at various times. We Caneisce 
that the growth mode becomes very broad when the source term is strong. The 
PVPF moves quickly toward larger sizes after the source term vanishes. The 
Mumerical artifact discussed above 18s “expected to be weaker when a very 


narrow Gaussian source pulse is used to represent the particle formation 


mechanisms. 


Figure 4.8 presents the number concentration and volume fraction 
profiles corresponding to the two modes, as calculated using the bimodal 
integral solution. At the early stage (0 ¢ t << t,) N, accounts for almost all 
the total number concentration of the aerosol. At the region of maximum 
formation rate (t = t,) N, is appreciably higher than before but always 


smaller than N,. Both N, and N, peak near t,, with N, decaying very 
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rapidly thereafter. It is apparent that N, = N, when the source term has a 
maximum, Correepo tie to quasi equilibrium between the coagulation and 
particle inception eee €xpressead wy dN/dt = (. Comparing the 
volume fraction profiles for the two modes we conclude that even at the 
moment of peak formation rate, where $, achieves a maximum, $, accounts for 
almost all the total volume fraction, ie. $, = 9. At later times when the 
source term becomes very weak or negligible, N, accounts for the biggest part 


of N since the source cannot further support the loss of particles from the 


inception mode due to intense coagulation. Figure 4.9 presents. the 
variation of the geometric standard deviation of the growth mode Og with 
time. As we see a strong particle formation rate implies a very wide 


distribution function due to explosively intense coagulation. After the source 
term weakens substantially the value of Ogg converges to a value Og. * 0.3. 
The bimodal solution finally reduces to the monomodal one after the source 


term has terminated and the inception mode vanishes. 
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4.5 _Bimodal solution predictions for the soot aerosol _in the coannular 
ethene/air diffusion flame. 


The bimodal solution is now used to model the soot formation 
processes in the axisymmetric ethene/air diffusion flame under investigation. 
The formation, growth and transport of soot particles along the path line 
merging with the annular region at a height z=15 mm (see Fig. 3.22) is 
studied in detail. As discussed previously during the experimental investigation 
of the soot field, these particles are first formed between the flame sheet and 
the annular region. The particle inception stage occurs in the lowest regions of 
fuemiaiie and its Characteristic time 4s O01 tne, order Of | ms. The newly 
formed particles follow the path lines leading to the interior of the flame and 
at the same time they are convected upwards to higher axial locations. During 
this transport the particles grow through surface growth reactions, while 
agglomeration results in the reduction of their number concentration. Finally 
the particles move into the upper oxidation region of the flame where they 


are ultimately depleted under the attack of the oxidative species. 


According to Figure 3.22 the path line which merges with the annular 
region at a height z=I15 mm, crosses the flame sheet at an axial location z=4 
mm. Since the exit plane of the fuel nozzle was chosen as the initial location 
(t=0) for the calculation of the particle residence times, this event corresponds 
to t=l16 ms. Combining the information provided by Figure 2.6A and B, 
reproduced from refs. 86 and 81 respectively, we can calculate the residence 
times t corresponding to various axial locations z, for the particles being 


transported along the path line merging with r, at z=15 mm. The values Geet 


= roo 


and z are given in the following table 


TABLE 4.4 


Residence Time Data®* for Particles Ultimately Transported 
on the Annular Region 


z (mm) a - 10 15 20 30 40 


t (ms) 16 20 28 34 39 47 54 
* Data from ref. 81. 


This table suggests that the particles are created over a time interval 
which starts at an instant t2?16 ms. These particles follow a path that merges 
with the annular region at a time t=34 ms, as they are transported upwards. 
This implies that the transport of these particles, from the region where they 
are formed (z=4 mm) to radial positions corresponding to the annular region 
(z=15 mm), does not last longer than 18 ms. During this transport the particles 
grow and their surface growth rate is expected to be lower than that 
calculated from the HES of the thermophoretic method for particles 
corresponding to r=r.. As soon as the particle path merges with the soot 
annulus, the surface growth processes can be described by the growth cara of 
particles transported along the annular region. These rates, which were 


derived by EM observations described previously, are given in Table 3.1. 


The value of the surface growth rate at a height z=10 mm for the 
path line under investigation, is expected®® to be lower than that of 9.4x1075 
g/cm? s which corresponds to the annular region. Assuming a value of 7.5x1075 


g/cm? s, a second order polynomial fit of the surface growth rate along this 
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particular path line gives f(t) for t?20.5 ms as 


ftict Avs Rt 4 Cts [g/em? s ; t in sec] 


where A=-0.443x1073, B=0.030 and C=-0.410. The above equation was used 
during the numerical simulation, to describe the specific growth rate derived 
using the thermophoretic sampling method, for residence times t?20.5 ms, and 
f=0 for t<20.5 ms. 

As we mentioned previously, the bimodal solution assumes particles of 
spherical shape. The calculated values of f(t) were obtained from measurements 
of the actual soot aggregates examining a series of electron micrographs. Let 
us consider an aggregate of volume V, which consists of Nh primary particles 
of uniform diameter D,,. If we neglect the bridging of the particles, the total 


surface of the aggregate is given by 


sB 2 
S ae mN,D, 


a 


Considering the equivalent sphere of volume V, we calculate its surface area 


as 


2 
| 2 
SF = nN , D, 


We then calculate the ratio 


(4.82) 


This ratio, which can be very different from unity, is a result of the 
nonspherical character of the soot aggregates. It is used in the numerical 


simulation as a multiplicative factor of the calculated values of f(t), in order 
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to compensate for the smaller area of the aggregated particles when they are 


considered to have a spherical shape. 


The variables of the numerical simulation are the following quantities: 


t T, describing the Gaussian source pulse used to simulate the particle 


Pas 
inception mechanisms along the path line ultimately merging with the soot 
annulus at z=15 mm, N, giving the total number of particles per unit volume 
generated by the particle inception mechanisms along the same path, Dg, and 
Og. describing the log-normal size distribution of the incipient soot particles, 
and finally fp defined by Eq. (4.82) which accounts for the nonsphericity of 


the soot aggregates. 


The data acquired through the optical method provide information on 
volume fraction, particle size and number concentration of the soot aerosol for 
particles corresponding to the annular region. Figure 2.7 reproduced from ref. 
86, which describes the results of optical observations from the soot field of 


this flame, was used to retrieve the following information 


TABLE 4.5 


Selected Optical Data for Particles in the Annular Region 
of the Ethene Diffusion Flame 


t Z (0) D N 


(ms) (mm) (cm ?/ cm? ) Gan) (cm73) 
20 Simer3 x1 On 20 erat 
28 10 ~ 1.8x10a¢ 73 1.8x101° 
34 1S5rhei7 x107% 110 herrero?? 
39 20 WE Sx h0e 135 8.7x10° 
47 30 PFs 0x 103% 150 8.6x10° 


54 AQMESTOx [0-° 140 1.3x102° 
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ameo9 valties’* listed’! One the= above’ table “were “taken™” using ~ the 
extrapolated? curves®tof "Figure “272! labeled’ 3.85 “ec/s’” The “reliability of “the 
optically calculated quantities decreases from 6 to N. These quantities were 
calculated using Mie light-scattering analysis for polydisperse spheres, with a 
self-preserving size distribution function and a complex refractive index 
m=1.57-0.561. As we see this series of assumptions causes some questions related 


to the reliability of the data acquired through the optical method. 


The particle volume fraction profile $(t) is used as the most reliable 
quantity in our numerical simulations. The values shown in table 4.5 
correspond to the annular region of the flame at all heights. Thus they do not 
describe the particles whose path ultimately merges with Be for heights below 
z=15 mm. At a height de LOCH the peak value of ® was found to be 1.8x107° 
cm?/cm*. Using Figure 3.22 we estimate that the particles under investigation 
cross the axial location z=10 mm at a radial position 0.5 mm outside pa ee 
corresponds to a reduction of the related value of 9$, suggesting that at z=10 
mm (t=28 ms) the numerically calculated $ should therefore be lower than the 


One corresponding to the annular region. 


An additional constraint suggested by the optical measurements was 
imposed on the numerical simulation. Equation (2.1) was used to calculate Q,, 
as a function of residence time for the particles being transported on the path 
under investigation. These values of Q,, are valid only when the particles are 
in the Rayleigh regime. This means that the calculated values of Q,, have to 
be consistent with the optically determined ones for heights z<10 mm (t¢28 


ms), where the Rayleigh theory can describe the soot aerosol reliably. From the 
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optical data we obtain that the lowest Q,, signal that can be detected is of 
the order of 1x107? cm@str7!. Using Figure 3.22 and refs. 85 and 86 for the 
particles whose path merges with r, at z=15 mm, we obtain Qericl0¥ cm7!str7} 
and Q..53-t) Og: cm7!str-! corresponding to heights z=5 mm and z=10 mm 
respectively. The numerical simulations of the soot aerosol using the bimodal 
integral solution have to be consistent with these two conditions in order to be 


acceptable. 


A Gaussian source pulse creating almost monodisperse particles 
(dg, =0.1) is used to simulate the particle inception mechanisms along the path 
ultimately merging with the soot annulus. The incipient soot particles are small 
(Dg, =1.5 nm), and are generated over a short time interval. The time of 
maximum naericl formation rate (t_) was chosen to correspond to a position 
close to the flame sheet on the fuel rich region. The total number of particles 


per unit volume generated by the source pulse was of the order of 10 Seems. 


We assumed a soot density of 1.86 g/cm? 


according to refs. 19,64,67 and a 
constant temperature T=2000 K for the soot aerosol. The value of f, was 
found to be 3.7 in order to account for the sharp increase of $ over the range 
10 to 30 mm of heights above the burner mouth. The following set of 
variables was used as one that gives agreement between the 9 profiles 


calculated using the bimodal solution, and the values determined by the optical 


method: 


t =17.5 ms, T,=0.6 ms, N,#1.2x10" cm™, D, =1.5 nm, 


og =0.1, and f,=3.7 


Bo 


The above set of variables satisfied the constraint of compatibility of the 


calculated Q,, values with the ones determined optically. 
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The assumption of no growth on the particle inception mode is now 
checked in order to assure the validity of the simulation. The bimodal 
character of the solution is preserved until the time t=23.5 ms, where the 
transition to the monomodal configuration occurs. At this instant it is 
N,/N, =500. The surface growth rate achieves non zero eaitiee after t=20.5 ms 
when N,/N,#=22. This means that when the surface growth term enters the 
calculation, the bimodal distribution is such that less than five percent of the 
particle population belongs to the inception mode. At times t=21 and 22 ms the 
ratio N,/N, achieves values 36 and 90 respectively. The rapidly decreasing 
ratio N,/N, during the time period 20.5<t<23.5 ms, illustrates the insignificant 
error introduced when we neglect the surface growth mechanisms acting on the 


particle inception mode. 


Figure 4.10 shows the particle volume fraction profile calculated by 
‘the bimodal solution which was used ‘to fit the profile determined 
experimentally using the optical tests. The particle generation mechanisms 
create over a short period of time an initial soot loading of about 
2.x10°8cm?/cm?, which is very low eeciouead to the loading higher up in the 
flame. Surface growth reactions cause a significant increase of the soot 
volume fraction thereafter, up to the higher portions of the flame where the 
reactivity of the particles decreases drastically and these reactions terminate. 
This is the point where oxidation mechanisms start to dominate thus causing a 
gradual decrease of $ with time. The oxidation mechanisms have not been 


modeled in the bimodal solution. 


Figure 4.11A shows the calculated particle diameter D,, versus 


residence time. The values found using the optical and  thermophoretic 
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techniques are also shown. We can see that the predictions of the model are 
consistently lower than the experimental values. Figure 4.11B shows the 
particle number concentration, as calculated by the bimodal solution, to achieve 
a maximum at tst,. The number concentration subsequently decreases 
monotonically with time as a result of meriithies mechanisms. The optical 
data corresponding to intermediate heights of the flame are also displayed. 
Figure 4.12A displays the variation of fy as calculated by the bimodal 
solution. The dimensionless second moment achieves a maximum at the time of 
maximum particle inception rate (t_), and subsequently decreases because of 
the intense surface growth mechanisms. As the surface growth rate weakens at 
later times, fy, converges slowly to the value corresponding to the self 
preserving distribution function. This value cannot be achieved in the short 
time available because of the presence of particle inception or surface growth 
processes throughout the lower and intermediate regions of the ethene diffusion 
flame. Figure 4.12A also demonstrates the numerical result that higher values 
of f, are associated with intense particle generation mechanisms. Low values 
of f, result due to surface growth reactions which tend to Bee: the size 
distribution. This is illustrated at the time interval during which the surface 
growth mechanisms are active, thus causing the values of f, to drop below the 
value corresponding to the self-preserving form. The coagulative mechanism of 
growth tends to broaden the distribution causing an increase of fy. The 
simultaneous presence of both growth processes creates a competition that 
results in the dominance of the surface growth reactions. However it is clear 
that the values of f, throughout the period of intense surface growth reactions 


are not very different from the self-preserving value. The result of “the 


above numerical simulation is in agreement with the studies of Bockhorn et al? 
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on low pressure hydrocarbon/oxygen premixed flat flames. In these flames they 
estimated the duration of particle inception to be 2 to 3 ms assuming an 


incipient particle diameter of 1.5 nm. 


Figure 4.12B shows the particle inception pulse S, used during the 
numerical simulation. We see that a maximum inception rate of about 1.1x10%6 
cm7%s"! was employed in order to simulate the first stage of the soot related 
processes in the ethene diffusion flame. The e~! half width of the pulse 
T,=0.6 ms and its value affects the maximum value of S, A higher value of T, 
would decrease the corresponding values of S,. A favorable reproduction of the 
experimentally determined volume fraction profile, can only be obtained using 


particle inception rates above the value of 10!® cm7* s-!. A sensitivity analysis 


related to the values of the particle inception rate was also performed. Using 


an incipient particle size of D, =1 nm we obtained an equally good fit of: the 
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“3 and the rest of the variables 


volume fraction profile with N,=4.2x10!% cm 
unchanged. During this simulation a maximum particle inception rate of 
3.9x101§ cm? s7! was employed. Once again the numerically calculated values 


of Q,, were consistent with the optically determined ones. 


Figure 4.13 shows the Q,, time variation corresponding to the lower 
portion of the flame as calculated by the bimodal solution using Rayleigh 


scattering theory. As we see the Q,., signal would not be detectable for t<20 


ms. This means that the optical method can not "see" the particle inception 
stage but starts detecting the particles at a later stage when they have grown 
over a certain size. This numerical result demonstrates that the particle 
concentration increase in the region of the flame sheet can be associated with 


the formation of large numbers of small particles (incipient particles). The 
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existence of the particle inception pulse, that is inferred by the application of 
the aerosol dynamic equation as depicted in Figure 4.12B, is seemingly 
impossible to verify experimentally. At this height, z=3-5 mm, the optical 
experiment does not yield results because the transmission is too low to permit 
any measurement in a reliable manner. The thermophoretic experiment lacks 
the spatial resolution that is necessary to identify a particle inception pulse, 
which is marked by coexisting coagulative growth at early times and surface 
growth soon thereafter. Thus the particle inception pulse is inferred largely by 
a process of numerical extrapolation’ at earlier times. While this procedure 
brings into question the accuracy of the derived inception rate, the importance 


of these attempts to understand the early events cannot be doubted. 


Our studies wml that the soot inception rate for the ethene 
diffusion flame is almost three orders of magnitude higher than that 
determined by Wersborg et al.°® for a low pressure acetylene/oxygen premixed 
flat flame. In that flame they observed a maximum soot volume fraction 
smaller by a factor of 500 (=1.5x107°? cm?/cm*) while the incipient particle 
size was assumed to be 1.5 nm. The maximum number concentrations were of 
the order of 101° cm-*. Bockhorn et al.? reported a particle inception rate of 
the order of 1014 cm-* s-! for their low pressure premixed flat flames 
assuming an incipient particle diameter of 1.5 nm. Harris et al.4? calculated an 


1016 cm? st for an atmospheric pressure premixed 


inception rate of 
ethene/argon/oxygen flat flame. They considered a particle inception pulse with 
a duration of 2-3 ms generating particles of the same size (1.5 nm), and they 
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observed number concentrations as high as 10/2 cm In their calculation 


they used a constant coagulation coefficient as Wersborg et al. did. Kent et 
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al5! determined an inception rate of 1015 cm-* s-! for an atmospheric pressure 
laminar ethene diffusion flame on a Wolfhard-Parker burner. As we see the 
inception rates calculated through our numerical simulations are generally in 
agreement with the rates reported previously for the ethene flames. The variety 
of the Bare conditions, not to mention the uncertainties involved in the 
calculations, could certainly account for the different values of particle 


inception rate that are reported. 


As described previously the quantity Re is defined as the ratio of the 
integrated particle volume fraction due to surface growth, over the integrated 
volume fraction due to particle formation. This quantity expresses the relative 
contribution of the two growth mechanisms to the final soot loading inside the 
ethene diffusion flame. The numerical simulation gave a value of R,=500 
corresponding to a height z=40 mm (onset of oxidation). This means that only 
0.2% of the final soot volume fraction is due to the particle inception 
mechanisms, while 99.8% is due to soot surface growth reactions. This result is 
in agreement with Sete ee performed by Harris and Weiner®®*9-4! 
(atmospheric pressure premixed ethene-air flat flames), and Bockhorn et al.? 
(low pressure premixed hydrocarbon-oxygen flat flames) who concluded that 


surface growth accounts for the bulk of soot generated in the flames they 


studied. 


Another significant result of the numerical simulation is the value of 
fp which was used as a factor to account for the nonsphericity of the soot 
aggregates and was defined by Eq. (4.82). A value of fp=3.7 was used 


throughout the flame. According to Eq. (4.82) this suggests an average number 


of primary particles per aggregate N,290. The average value determined by the 
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thermophoretic method for particles corresponding to the annular region at a 
height z=15 mm, was Nase which is in fair agreement with the value 
obtained numerically. However this expresses an average value and is not 
expected to describe accurately all the local values of Né inside the flame. 


Howard, Wersborg and Williams*® 


estimated that N,=10 for a low pressure 
premixed acetylene/oxygen flat flame. The higher number concentrations in our 
diffusion flame suggest higher agglomeration rates which imply a larger value 


of N,. 


Equation (4.34) of the bimodal solution was employed to calculate the 
surface area of the soot particles ultimately transported on r, as a function of 
residence time. As shown in Figure 4.14A, during the time interval of intense 
particle inception, an increase of the surface area occurs which is related to 
the generation of a large population of fine particles. The soot surface area 
subsequently decreases because of strong coalescence between the particles. This 
decrease occurs only in the absence of surface growth which tends to increase 
surface area. After the SE of the surface growth mechanisms (the specific 
growth rate used is shown in Figure 4.14B) there is a strong competition of 
the two growth processes. The steady increase of the surface area, displayed by 
Figure 4.14A during that period, suggests the dominance of the surface growth 
term over the coagulation term. The subsequent cessation of the surface growth 
mechanisms triggers another stage of surface area decrease, associated with 
slower coagulation mechanisms as a result of the lower particle number 
concentrations. The above picture is not consistent with the observations of 
Harris and Weiner?®°%4! for a premixed ethene/air flame under atmospheric 
pressure conditions. They reported an almost constant value of the soot surface 


area for three different stoichiometries. A series of assumptions employed by 
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both methods does not permit any conclusions about the reliability of the data, 


which refer to two flames of different nature. 


The number concentration increase on either side of r., as determined 
by the optical method, was also investigated using the bimodal solution. As we 
mentioned previously, the annular region coincides with the location where the 
surface growth mechanisms display their maximum intensity. The surface 
growth reactions on both sides of the annular region (r>r, and r<r,) can be 
substantially weaker as demonstrated by the acetylene production rate studies 


performed by Miller®. 


This behavior was modeled using the bimodal solution 
with no surface growth term. The rest of the variables were the ones used 
previously to describe the soot aerosol on the path line merging with r, at 
z=15 mm. An important point was inferred from these calculations. Figure 4.15 
shows the time variation of ‘the particle number concentration and the particle 
size with and without the surfacé growth term. As we see the absence of 
surface growth reactions results in slower coagulation rates leading to higher 
number concentrations and smaller particle diameters. This numerical result 
illustrates that the number Shitcentrahes anomaly, as determined by the optical 


measurements, may be a result of the weaker surface growth reactions occuring 


on both sides of the annular region. 
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CHAPTER 5 


CONCLUSIONS AND RECOMMENDATIONS 


A technique for extracting soot or other combustion generated (fumed) 
particulates from flames has been devised that affords a detailed examination 
of particle morphology and its regional variation within a flame environment. 
The technique is useful for gaining information on particle morphology that is 
not yielded by scattering and extinction measurements. The thermophoretic 

2 


technique provides additional information for absorbing particles (Q,, = nv’, 


K.., * nv), while for dielectric particles (Q,, & K,, * nv’) the method is 


ext 
even more important. The method uses no assumption on the shape or the size 
distribution of the particles examined and provides additional information on 


the particle field which is necessary to reduce the data yielded by the optical 


experiments. 


The method was used to study an overventilated coannular laminar 
ethene/air diffusion flame that releases no soot to the surroundings. The 
morphological features of the particles, examined by TEM for best image 
quality, provide not only valuable qualitative information on _ particle 
agglomeration, surface growth, and oxidation but also quantitative data on 
primary particle size as a function of flame coordinates. The technique proved 
to be rather unsuccessful in the lowest portion of the flame where particle 
inception occurs and where the reproducibility of the morphological data 
showed no consistency. Observations of soot samples collected from the low and 
intermediate portions of the flame consistently show an increasing primary size 
with height above the burner. This proves the dominance of surface growth 


reactions on the morphology of soot particles in the ethene flame 
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environment. The specific surface growth rate of particles transported on the 
annular region can be calculated throughout the intermediate regions of the 
flame without any approximations on the form of the particle size distribution, 
if particle velocity data are also available. A maximum growth rate of 
10.5x10°> g/cm?s was measured occuring at an axial location z=15 mm above 
the burner mouth. The primary particle size was found to decrease as the 
particles move to the upper regions of the flame, probably through oxidative 


attack of various species. 


Our studies further show that, at a Specific axial location for 
intermediate heights in the flame, the primary particle size and the degree of 
agglomeration peak on the radius of maximum soot volume fraction (r=r,). 
Both SEM and TEM observations show that the primary particle size and the 
state of agglomeration are substantially reduced on the air side of r., where 
particles smaller than 5 nm were detected throughout the interval 10<z<40 mm. 
A similar trend is also present on the fuel side of r, but the particle sizes 
there are not as small. However the observations corresponding to r<r, are less 
reliable since these regions are surrounded by an Ren is of high soot 
concentration. The above observations led to the conclusion that the soot 
particles are formed locally in low regions of the flame near the reaction zone 
(maximum temperature region). This mechanism can account for the increase 
of the particle number concentration close to the reaction zone as determined 


by the optical measurements. 


A comparison of the volume equivalent mean diameters obtained by 
the thermophoretic and the optical methods at a height z=l15 mm above the 


burner mouth showed a fair agreement between the two techniques, thus 


as 


supporting the assumptions adopted by the optical tests for the data reduction 
at this height. The number of soot aggregates manually analyzed using an 
image analysis system proved to be inadequate to determine the form of the 


size distribution in a reliable manner. 


The bimodal integral solution, used to model the soot formation 
processes in the ethene flame, provided additional information on the transpozt 
of particles along the path line merging with the annular region at a height 
z=15 mm. Our numerical simulations assumed incipient particles of 1.5 nm 
diameter generated by a Gaussian source pulse with ae! half width of 0.6 
ms, creating a soot loading of the order of 107° cm?/cm’,. A maximum 


particle inception rate of 101° cm? s-! 


was employed, a rate that is in 
agreement with the values reported previously for ethylene flames. The amount 
of soot generated during the particle inception stage was found to be 
negligible (0.2%) compared to that caused by surface growth reactions which 
accounts almost entirely for the bulk of soot ultimately formed. The value of 
the dimensionless second moment fy was found to be higher during the 
particle inception stage and was almost constant during the interval of intense 
surface growth. Finally the soot surface area was found to increase during the 


surface growth stage which lasts throughout the low and intermediate regions 


of the flame. 


Some recommendations on future directions of the current program 
are now suggested: 
The investigation of the higher portions of the flame can provide 


valuable information on the soot oxidation rate in the ethene flame. The 
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potential of the method to obtain surface growth or oxidation rates has been 
proven and can be used to study the phenomena that trigger the release of 


soot to the environment under appropriate flow conditions. 


Even though the soot field of the coannular ethene diffusion flame 
has been thoroughly investigated, though not completely understood, no 
extensive study of the chemical structure of the flame is available. Some 
previous attempts have given only limited amount of data proving the 
difficulty of such measurements. A mass spectrometric analysis of the basic 
species profiles would provide valuable information about the chemical steps 


leading to the formation, growth and depletion of soot throughout the flame. 


In the theoretical model some improvements can be incorporated in 
the formulation of the collision term of the dynamic equation in order to 
‘account for the more complicated agglomeration mechanisms. An appropriate 
collision factor, modeling the nonspherical character of the soot aggregates, 
would obviously create a better description of the collision mechanisms of the 


particles inside the flame. 


The thermophoretic sampling method can be used to provide soot 
morphological data related to other diffusion or premixed flames that have 
been already investigated optically and chemically. The combination of these 
techniques can contribute potentially to a detailed understanding of the soot 


formation processes in flames, which is far from our grasp at this time. 
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APPENDIX _A 


Numerical evaluation of the collision integrals 


The collision integrals were defined by Eqs. (4.7). Assuming a_ log-normal 


form of the PVDF and using the following transformation 


In V/Ve. 

xX, = o Oen ae (hawt MP4) 
= 
BYa Og. 
In u/v, 

X, = — ee (i) -wenle 2) 
= 
3°%2 Og. 

it can be shown that 
(0) 
Z i 
fcivnas G. Abies 4 exp (-x,. Xo "dx, dx, (A.1) 


and 


g. 
Z¢m, = pla [ f, exe Xo ]G, le ines ee exp (-x, - x, dx, dx, (A.2) 


where pin xh = exp| 32 oa. (x, + 7] 
1 
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OF se : 
In order to express Zfm, in the same form, define 


y | ln w/v 
€ Mie = 
Ve 
3 oe aey 
In ihe 4s 
Ye 
= 
3M) oF 
Then 
i i o @ 
Die ee 21 2 2 
“fm, = 5 aa satin al Gye Yq) exp (-y, - Yo }dy, dy, (A.3) 
(2) ate ae: 
Aon oad Bi ‘Te | fy (ways) G, (v1. ¥3] exp (-y, = J }dy, dy, (A.4) 
where 
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As we see from (A.1) - (A.4) all collision integrals Zp, ; 1 = 1, 2) Dogan 
1 


k = 0, 2 are 61 thepeeneral { orm 


eo 


(k) 
Zem, = F(x, y) exp (-x? - y? dx dy (A.5) 


8 
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This form of integral can be readily evaluated using the Gauss-Hermite 


Quadrature (GHQ) rule!?5 expressed by 


ed n 
eX” f(x) dx = ) vaeibeat (A.6) 
= m=] 


where w_: weight factors 


xX, zeros of Hermite polynomials 


The values of w_, x_ @orizg < n < 20 ane given in Table 2:4 oterei ai 


For high accuracy we chosetn = 20. 
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Applying the GHQ rule on a double integral we can express (A.5) in a Series 


form given by 


Nn n 
(k) 
= =] 


m=1] j 


(A.7) 


: k 
This formula was used to evaluate numerically the integrals ari {rare D 
1 


aie Ko) 0, 72. 


(k) 


The values of Z (t= 12 and kK = 0.2) "ican ben closely “approximated by 
fm, 


interpolation formulae. Dobbins and Mulholland?* found that 
k 
OWS) o,, < 1.2, Zee can be approximated by 
1 1 
(0) ie Pini er 
Zz | ZY ex _ fo} 
fm; "8 on aes 
(2) le Stee 
ss v> 
Tinh Seems exe (Fp ey, | 


However the accuracy is not acceptable for values of a, > 1.2. 


1 
(k) 
values of 0, we approximate Zp, by 
1 1 


1 
(0) e 6 2 i) 
=- a.+ a,o tag O. eee te st ohh 2 
fm, av 2 eral SLs Lov. ek 
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in the range 


(A.8) 


(A.9) 


For these 


(A.10) 


(A.11) 


where 


rs) 
il 


~8.346 ox 10%, pa, 91.842 x 107) as an 207 ome 


and 


ra) 
rT] 


, = -2.255 x 10, a, = 1.359) agenn6.306x 10° 


These interpolation formulae provide good accuracy for high values of o, 
1 


(see table A.1) 
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TABLES AL 


Comparison of Collision Integrals with Interpolation Formulae 


oy, oy - a GHQ* Eq.(A.10)  GHQ Eq. (A.11) 
1.50 0.50 550 1.557 21.536 21.461 
1.80 0.60 1.941 1.944 85.063 85.242 
2.10 0.70 2.568 2.574 440.44 442.08 
2.25 Ov5e - : 1,112.5 1,115.4 

2.40 0.80 3.645 3.646 3,016.8 3,017.5 
2.55 0.85 4.471 4.467 8,785.8 8,763.5 
2.70 0.90 5.599 5.588 Zr 21357 
2.85 0.95 7.161 7.151 92,199 91,920 
3.00 1.00 9.354 9.375 331,307 332,898 


* Evaluated using the Gauss-Hermite Quadrature rule expressed by Eq. (A.7) 


ot Es 


REFERENCES 


1. ASTM:D3849-85. Standard test method for carbon black-primary aggregate 
dimensions from electron microscope image analysis. Annual Book of ASTM 


Standards. 


2. Baumgartner, L., Hesse, D., Jander, H., and Wagner, H. Gg. (1984). Rate 
of soot growth in atmospheric premixed laminar flames. Twentieth Symposium 


(International) on Combustion, The Combustion Institute, pp.959-967. 


3. Bockhorn, H., Fetting, F., and MHeddrich, A. (1986). Investigation of 
particle inception in sooting premixed hydrocarbon-oxygen low pressure flames. 
Twenty First Symposium (International) on Combustion, The Combustion 


Institute, in press. 


4. Bockhorn, H., Fetting, F., Heddrich, A., and Wannemacher, G. (1984). 
Investigation of the surface growth of soot in flat low pressure hydrocarbon 
oxygen flames. Twentieth Symposium (International) on Combustion, The 


Combustion Institute, pp.979-988. 


5. Bockhorn, -H:;~-Fetting;~F""Meyer, "U:,” Reck;” R.,” and” Wannémacher.ec:, 
(1981). Measurement of the soot concentration and soot particle sizes in 
propane oxygen flames. Eighteenth Symposium (International) on Combustion, 


The Combustion Institute, pp.1137-1147. 


6. Bockhorn, H., Fetting, F.. Wannemacher, G., and Wenz, H. W. (1982). 
Optical studies of soot particle growth in hydrocarbon oxygen flames. 


Nineteenth Symposium (International) on Combustion, The Combustion 


= 30% 


Institute, pp.1413-1420. 


7. Bockhorn, H., Fetting, F.. and Wenz, H. W. (1983). Investigation of the 
formation of high molecular hydrocarbons and soot in premixed hydrocarbon- 


oxygen flames. Ber. Bunsenges. Phys. Chem. 87,1067. 


8. Boedeker, L. R., and Dobbs, G. M. (1986). - CARS temperature 
measurements in annular soot zones of axisymmetric laminar diffusion flames. 
Paper 36, Technical meeting, Eastern Section of the Combustion Institute, San 


Juan, Puerto Rico. 


9. Burke, S. Pang” Schumann. o.- Ey W. = (1928). Diffusion 


flames. Industrial Engineering Chemistry 20,998. 


10. Calcote, H. FF. (1981). Mechanisms of soot nucleation in flames-A 


critical review. Combustion and Flame 42,215. 


ieeeecoarnahnan. 6. Luther, H.- A:; and Wilkes, J.. ©O. (1969). Applied 


Numerical Methods, John Wiley & Sons, New York. ~ 


12. Chippett, S. and Gray, W. A. (1978). The size and optical properties of 


soot particles. Combustion and Flame 31,149. 


PemeeGoées DD. PSceHaynese:B. WSusandesteinteldael nl(1981). eldentificationsof 
a source of argon-ion-laser excited fluorescence in sooting flames. Combustion 


and Flame 43,211. 


14. D’Alessio, A. (1981). Laser light scattering and fluorescence diagnostics 
of rich flames produced by gaseous and liquid fuels. Particulate Carbon 
Formation during Combustion, D. C. Siegla and G. W. Smith (Eds.), Plenum 


Press, New York, p.207. 


-139- 


15. D’Alessio, A., Beretta, F., and Venitozzi, C. (1972). Optical investigations 


on soot forming methane-oxygen flames. Combustion Sci. and Technology 5,263. 


16. D’Alessio, A., DiLorenzo, A., Borghese, A., Beretta, FE. andeeNigeuee 
(1977). Study of the soot nucleation zone of rich methane-oxygen flames. 
Sixteenth Symposium (International) on Combustion, The Combustion Institute, 


pp. 695-708. 


17. D’Alessio, A., DiLorenzo, A., Sarofim, A. F., Beretta, F., Masi, S., and 
Venitozzi, C. (1975). Soot formation in methane-oxygen flames. Fifteenth 
Symposium (International) on Combustion, The Combustion Institute, 


pp.1427-1438. 


18. Dannenberg, E. M. (1952). Carbon black dispersion and reinforcement. 


Industrial and Eng. Chemistry 44,813. 


19. Dannenberg, E.. M.,. Heckman, F., .A., and Medalia, .A.s. ego) 
Structure and particle size of carbon black. Presented at the International 


Rubber Conference, Paris, France; Cabot Corporation research paper 57-88. 


20. Dave, J. V. (1978). Subroutines for computing the parameter of the 
electromagnetic radiation scattered by a sphere. Report 320-3237, IBM, Palo 


Alto Scientific Center, Palo Alto CA. 


21% Dobbins, R. A., and Megaridis, C. M. (1986). Soot morphological 
characteristics by means of thermophoretic sampling. Poster paper PS-009, 
Twenty First Symposium (International) on Combustion, The Combustion 


Institute, Munich, W. Germany. 


-140- 


eZee WoovInS,. R. -A;) and .Megaridis, .C. M..(1987). Morphology of 
flame-generated soot as determined by thermophoretic sampling. Langmuir 3, 


No. 2, 254. 


23. Dobbins, R. A., and Mulholland, G. W. (1984). Interpretation of 
optical measurements of flame generated particles. Combustion Sci. and 


Technology 40,175. 


24. Eisner, A. D., and Rosner, D. E. (1985). Experimental studies of soot 
particle thermophoresis in nonisothermal combustion gases using thermocouple 


response techniques. Combustion and Flame 61,153. 


25. Ferziger, J. H. (1981). Numerical Methods for Engineering Application, 


John Wiley & Sons, New York. 


#45 Flower, W. bE, and “Bowman, —C. T. (1984). Measurements of the 
structure of sooting laminar diffusion flames at elevated pressures. Twentieth 
Symposium (International) on Combustion, The Combustion Institute, 


pp.1035-1044. 


eye) Flower, WwW." Ly and Bowman, Cr" TT.” (1986)> Soot production in 
axisymmetric laminar diffusion flames at pressures from one to _ ten 
atmospheres. Twenty First Symposium (International) on Combustion, The 


Combustion Institute, in press. 


28. Forsythe, G. E., Malcolm, M. A., and Moler, C. B. (1977). Computer 


Methods for Mathematical Computations, Prentice Hall, Englewood N.J. 


29. Frenklach, M., and Harris, S. J. (1987). Aerosol dynamics modeling 


using the method of moments. J. Colloid and Interface Sci., in press. 


wt Sad 


30.° Friedlander; S.~ K. (1977). Smoke, Dist and Haze, John” Wiley sommoaum 


New York. 


£08 Friedlander, “S: Ke, sand 8Wage ec: S.) (1966). The _ self-preserving 
particle size distribution for coagulation by Brownian motion. J. Colloid and 


Interface Sci. 22,126. 


32. Gelbard, F. (1984). MAEROS: Computer program and algorithm review. 


Aerosol Sci. and Technology 3,117. 


33. Gelbard, F., and Seinfeld, J. H. (1980). Simulation of multicomponent 


aerosol dynamics. J. Colloid and Interface Sci. 78, 485. 


34. Glassman, I., and Yaccarino, P. (1981). The temperature effect in 
sooting diffusion flames. Eighteenth Symposium (International) on Combustion, 


The Combustion Institute, pp.1175-1183. 


35. Goren, S. L. (1977). Thermophoresis of aerosol particles in the laminar 


boundary layer of a flat platé. J. Colloid and Interface Sci. 61,77. 


36. Graham, S. C.,,and Homer, J. _ B. (1973). . Coagulation of .moltentdeaa 


aerosols. Faraday Symposium Chem. Soc. 7,85. 


37. Graham, S. C., and Robinson, A. (1976). A comparison of numerical 
solutions to the self-preserving size distribution for aerosol coagulation in the 


free-molecule regime. J. Aerosol Sci. 7,261. 


38. Harris, S. J., and Weiner, A. M. (1983). Surface growth of soot 
particles in premixed ethylene/air flames. Combustion Sci. and Technology 


BH yd isge', 


=[22= 


evemeriairis, Ss. J. and ‘Weiner, A. M.(J983))" Determination of the, rate 


constant for soot surface growth. Combustion Sci. and Technology 32,267. 


40. Harris, S. J., and Weiner, A. M. (1984). Soot particle growth in 


premixed toluene/ethylene flames. Combustion Sci. and Technology 36,10: 


41. Harris, S. J.. and Weiner, A. M. (1985). Chemical kinetics of soot 


particle growth. Ann. Rev. Phys. Chem. 36,31. 


a2 ueriarris; S/7)J Weiner.) AY *M.) Tand “Ashcraft, C. “C. (1986). Soot particle 


inception kinetics in a premixed ethylene flame. Combustion and Flame 64,65. 


a3-2>raynes, B. “S.. *Jander, *H:,) and” Wagner, H.’ Gg. (1980). Optical studies 
of soot-formation processes in premixed flames. Ber. Bunsenges. Phys. Chem. 


84,585. 


44>" Haynes, _B>~ SS: and: Wagner,* H. Gg: (1980). “ Sooting structure in a 


laminar diffusion flame. Ber. Bunsenges. Phys. Chem. 84,499. 


ase Haynesie RB. eSie and #Wagner, Ha" Ger (i9si)ai Soot formation! jocProg. 


Energy Combust. Sci: 7,229. 


46)"%Hess, WP M3) Ban} “Le” 'L:3’and):McDonald,"G: ¢."(1969). © Carbon* black 
morphology: I. Particle microstructure. II]. Automated EM analysis of aggregate 


size and shape. Rubber Chem. and Technology 42,1209. 


47. Homann, K. H., and Wagner, H. Gg. (1967). Some new aspects of the 
mechanism of carbon formation in premixed flames. Eleventh Symposium 


(International) on Combustion, The Combustion Institute, pp.371-379. 


—Va3- 


48. Howard, J. B., Wersborg, B. L., and. Williams. 3G, C+ (19733. 
- Coagulation of carbon particles in premixed flames. Faraday Symposium Chem. 


Soc. 7.109: 


49. Jacoda, I. Ji, Prado, Gowand Canaye, Ji(19s0). An experimental 
investigation into soot formation and distribution in polymer diffusion flames. 


Combustion and Flame 37,261. 


50. JANAF Thermochemical Tables (1970). U. S. Government Printing Office, 


Sec. Edition, NSRDS-NBS37. 


51. Kent, J. H., Jander, H., and Wagner, H. Gg. (1981). Soot formation in 
a laminar diffusion flame. Eighteenth Symposium (International) on 


_ Combustion, The Combustion Institute, pp.1117-1126. 


2. Rents H., and Wagner, H. Gg. (1982). Soot measurements in 


laminar ethylene diffusion flames. Combustion and Flame 47,53. 


53. Kent, J.. H.,.and Wagner,- H. «-Gg..(1984). . Why» do ediffusiony flames 


emit smoke? Combustion Sci. and Technology 41,245. 


54. Kunugi, M., and Jinno, H. (1967). Determination of size and concentration 
of soot particles in diffusion flames by a light-scattering technique. Eleventh 
Symposium (International) on Combustion, The Combustion Institute, 


pp.257-266. 


55. Lahave. ais, and, «Prados G. (1978). Mechanisms of carbon black 
formation. Chemistry and Physics of Carbon, P. L. Walker, Jr. and P. A. 


Thrower (Eds.), Marcel Dekker, New York, vol.14, p.167. 


-144- 


56. Lahaye, J., Prado, G., and Donnet, J. B. (1974). Nucleation and growth 
of carbon black particles during thermal decomposition of benzene. Carbon 


eae p 


Sampaio i5.,- Briedlander, S. K., Pich, J., and MidvaG. eMIK1972) The 
self-preserving particle size distribution for Brownian coagulation in the free 


molecular regime. J. Colloid and Interface Sci. 39,395. 


pompecce-hoe W. (1983). Change of particle size distribution during Brownian 


coagulation. J. Colloid and Interface Sci. 92,315. 


By. wlLieeekK. 1.5 oande! ChengH.a(1984)ig0 Coagulation @raté of) polydisperse 


particles. Aerosol Sci. and Technology 3,327. 


Gi Leetag KO Wea Chen, iH, 4 and 4Giesekesodiia7 Avs £1984). Log-normally 
preserving size distribution for Brownian coagulation in the free-molecule 


‘regime. Aerosol Sci. and Technology 3,53. 


61. Medalia, A. I. (1967). Morphology of aggregates- I. Calculation of shape 
and bulkiness factors. Application to computer-simulated random flocs. Sf 


Colloid and Interface Sci. 24,393. 


62. Medalia, A. I. (1970). Morphology of aggregates- VI. Effective volume 
of aggregates of carbon black from electron microscopy; Application to vehicle 
absorption and to die swell of filled rubber. J. Colloid and Interface Sci. 32, 


INoOeise1 15, 


63. Medalia, A. I. and Heckman, F.A. (1971). Morphology of aggregates- VII. 
Comparison chart method for electron microscopic determination of carbon 


black aggregate morphology. J. Colloid and Interface SEm OOlWING LD: 


=145- 


64. Medalia, A. I., Dannenberg, E. M, Heckman,.F. A., and Cotten; Giak: 
(1973). Characterization of new technology carbon blacks. Rubber Chem. and 


Technology 46, No.5, 1239. 


65. Medalia, A. I. and Heckman, F. At (1969). Morphology of aggregates- 
II. Size and shape factors of carbon black aggregates from electron microscopy. 


Carbon 7,567. 


66. Medalia, A. I., Heckman, FF. A., and, Harling ab: F. (1969). 
Morphology of aggregates- IV. Particle size and structure of carbon black 


from electron micrographs. J. Rubber Res. Inst. Malaya 22,381. 


67. Medalia, A. I. and Richards, L. W. (1972). Tinting strength of carbon 


black. J. Colloid and Interface Sci. 40, No.2, 233. 


68. Miller, J. H. (1986). The chemical structure of an atmospheric pressure 
ethylene/air diffusion flame: Preliminary results. Paper 33, Technical meeting, 


Eastern Section of the Combustion Institute, San Juan, Puerto Rico. 


69. Miller, J. H, Mallard, W. G. and Smyth, K. C. (1982.5 
observation of laser-induced visible fluorescence in sooting diffusion flames. 


Combustion and Flame 47,205. 


70. Miller, J. H., Mallard, W. G. and Smyth, K. C. (1986). Chemical 
production rates of intermediate hydrocarbons in a methane/air diffusion 
flame. Twenty First Symposium (International) on Combustion, The Combustion 


Institute, in press. 


71. Millikan, R. C. (1962). Sizes, optical properties, and temperatures of 


soot particles. Temperature, Its Measurement and Control in Science and 


=Tn0> 


Industry, Vol.3, part 2, Reinhold, New York, 497. 


dies Mitchell, . R. Es) SarofinineraA: Bee and ee ClOnroure, (ec. A. (1980). 
Experimental and numerical investigation of confined laminar diffusion flames. 


Combustion and Flame 37,227. 


73. Mulholland, G. W., Mountain, R. D., and Baum, H. (1986). Simulation 
of aerosol agglomeration in the free molecular and continuum flow regimes. 


NBSIR 86-3342. 


74, Neoh, K. G., Howard, J. B., and Sarofim, A. F. (1984). Epiect 
of oxidation on the physical structure of soot. Twentieth Symposium 


(International) on Combustion, The Combustion Institute, pp.951-957. 


75. Prado, G., Garo, A., Ko, A., and Sarofim, A. (1984). Polycyclic aromatic 
hydrocarbons. Formation and destruction in a laminar diffusion flame. 
Twentieth Symposium (International) on Combustion, The Combustion Institute, 


pp.989-996. 


eGoeerrado, G. Jacoda, J.,~Neoh, KK. and Lahavyes J, (1981). A: study of soot 
formation in premixed propane/oxygen flames by in-situ optical techniques and 
sampling probes. Eighteenth Symposium (International) on Combustion, The 


Combustion Institute, pp.1127-1136. 


77. Prado, G. and Lahaye, J. (1981). Physical aspects of nucleation and 
growth of soot particles. Particulate Carbon Formation during Combustion, D. 


C. Siegla and G. W. Smith (Eds.), Plenum Press, New York, p.143. 


Teemerrano. Gy P+ Le¢. Me L.. nites, Roe AS Houlte Der and, Doward, J. B. 


(1977). Soot and hydrocarbon formation in a turbulent diffusion flame. 


-147- 


Sixteenth Symposium (International) on Combustion, The Combustion 


Institute, pp.649-661. 


79. Rosner, D. E. (1985). Mass transfer across combustion gas thermal 
boundary layers-Power production and materials processing implications. Invited 


lecture, 23rd ASME/AIChE National Heat Transfer Conference, Denver CO. 


80. Saito, K., Williams, F. A., and Gordon, A. S. (1986). A study of 
the two-color soot zone for small hydrocarbon diffusion flames. Combustion 


Sci. and Technology, in press. 
81. Santoro, R. J. (1984). Personal communication. 


82. Santoro, R. J., and Miller, J. H. (1987). Soot particle formation in 


laminar diffusion flames. Langmuir 3, No. 2, 244. 


83. Santoro, R. J., and Semerjian, H. G. (1984). Soot formation in 
diffusion flames: Flow rate, fuel species and temperature effects. Twentieth 
Symposium: (International) on Combustion, The Combustion Institute, pp. 


997-1006. 


84. Santoro, R. (ON ET Vel G., and Dobbins, R. A. (1983). 
Interpretation of optical measurements of soot in flames. AIAA Eighteenth 


Thermophysics Conference, Montreal, Canada. 


85. Santoro, R. J., Semerjian, H. G., and Dobbins, R. A. (1983). Soot 


particle measurements in diffusion flames. Combustion and Flame 51,203. 


86. Santoro, R. J., Yeh, T. T., and Semerjian, H. G (1985); The) tragsaon 
and growth of soot particles in laminar diffusion flames. 23rd ASME/AIChE 


National Heat Transfer Conference, Denver CO. 


a T= 


87. Seigneur, C., Hudischewskyj, A. B., Seinfeld, J. H., Whitby, K. T., Whitby, 
E. Baste DLOCK, Mem) Rip and “Barnes: ete NL (1960), ~=«Sitmulation. of aerosol 
dynamics: A comparative review of mathematical models. Aerosol Sci. and 


Technology 5,205. 


88. Skolnik, E. G., and McHale, E. T. -*(1980): Sampling of soot in 


diffusion flames. Combustion and flame 37,327. 


89. Smyth, K. Gyetand = Miller, J. He” «©6(19686); = Soot) “inception in 
hydrocarbon diffusion flames. Invited paper C, Technical Meeting, Eastern 


Section, The Combustion Institute, San Juan, Puerto Rico. 


SOSeromythekauG. Millers: J. H., Dorfman)’ R. “°C *Mallard, “Wo "’G..”~ and 
Santoro, R. J. (1985). Soot inception in a methane/air diffusion flame as 


characterized by detailed species profiles. Combustion and Flame 62,157. 


91. Svehla, R. A. (1962). Estimated viscosities and thermal conductivities 


of gases at high temperatures. NASA Technical Report R-132. 


po ewalbou Le Cheng.7 Rivi@Ke Schefers*R. ®W.Pand Willis, ~Di°R: “(1T980). 
Thermophoresis of particles in a heated boundary layer. J. of Fluid Mechanics 


LOT 737: 


93. Tesner;* P. A. (1959). Formation of dispersed carbon by thermal 
decomposition of hydrocarbons. Seventh Symposium (International) on 


Combustion, The Combustion Institute, pp.546-553. 


94. Tesner, P. A. (1973). Formation of soot particles. Faraday Symposium 


Chem. Soc. 7,104. 


=149= 


95., Tesner, P. .A., Snegiriova,:d; D.,,and,.Knorre;. Vii G, (19ddosepRimetice 


of dispersed carbon formation. Combustion and Flame 17,253. 


96. Wagner, H. Gg. (1981). Soot formation-An overview. Particulate 
Carbon Formation during Combustion, D. C. Siegla and G. W. Smith (Eds.), 


Plenum Press, New York, p.1. 


97. Warren, D. R., and Seinfeld, J. H. (1985). Simulation of aerosol size 
distribution in systems with simultaneous nucleation, condensation and 


coagulation. Aerosol Sci. and Technology 4,31. 


98. Wersborg, B. L., Howard, J. B., and Williams, G. C. (1973). Physical 
mechanisms in carbon formation in flames. Fourteenth Symposium 


(International) on Combustion, The Combustion Institute, pp. 929-940. 


99. Whitby, Ke T.2(198)5: Determination of aerosol growth rates in the 


atmosphere using lumped mode aerosol dynamics. J. Aerosol Sci. 12,173. 


100. Wolfhard, H. G., and Parker, W. G. (1949). A new technique for the 
spectroscopic examination of flames at normal pressures. Proc. Phys. Soc., 


London, Sect. A, 62,722. 


LOL Wolfhard, H. Gi and Parker. |W: G. ((1952); A spectroscopic 
investigation into the structure of diffusion flames. Proc. Phys. Soc., London, 


Séctu,AmGow, 


2750 


Luminous 
envelope 


Height 


Blue zone 


Fuel 


Figure 2.1 Typical structure of a nonsmoking coannular laminar diffusion 
flame. 
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Figure 2.2 (a) Soot volume fraction, (b) particle diameter, and (c) number 
concentration as a function of z and r for the ethene’ diffusion 


flame. Reproduced from ref. 85. 
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Figure 2.3 Velocity profiles in the ethene diffusion flame; (a) axial velocity 


component, (b) radial velocity component. Reproduced from ref. 86. 
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Figure 2.4 Temperature profiles in the ethene diffusion flame at various axial 


locations z (uncorrected thermocouple measurements). Reproduced from ref. 82. 


-154- 


‘poXvjdsip ose st wojsks 9}euIps009 


MYL “ULTJ UOISNJJIp ivurMe, svjnuuLos B& JO dIN}ONI}IS JUOIy OWT GZ INS 


yall 
| ; | 
Pe | 
| | 
| apis yo jens apis ude] jen4 
| 
| 


if 
/ 
(91) uolBeu sojnuuy / 


juol} SnoulWN > 


(9U0Z UO!}9DA1) 
JUS BUWD]4 


-155= 


[g pue 9g ‘SJol WoOly poonpoidoy 


AWOOJOA 94} WOJJ payepnojeos ‘sown 


(WI) “SOd WIOBY 
eK) S20 00°0 = S2"Q- OS 'Q- = SL “Qe 


(a) 


“AJOATIOIdSII 
‘QWIL]J UOISNyJIp dudyIs Oy} JOJ SJUNWOINSvOW 


asuapisol (gq) pue ‘syed ojonieg (VY) OZ D4INdILy 


(WW) NOILISOd 1IVIGVY 
Ol g 0 G- 


dy sauing 


00°t+, 


00°02 00° 


00°0h 


00°09 


00°08 
SWIL 33N301S3Y 


oa°00s 


ao°det 
(SW) 


00°095 00°Ons 


00°08s 


Oc 


OV 


O09 


08 


OOt 


(WW) JONVLSIO VIX 


=1>0= 


100 


1078 


ee 
O 
= 
O re) 10 
ps \ 0 20 40 60 80 =—:100 
a : TIME (ms) 
od 2.3 cc/e ° 
Z| 1078 \ 
> H 
Pa 13.66 cc/e 
Oo ' ae) 10'! 
O ; ' 
n \ E 
b — 
: z o-~*® 
= ° 74d cele 
re) = 
: 7) 
: z 
: WwW 
i = 10!0 
10a" - 
(0) 20 40 60 80 100 cD 
TIME (ms) 2 
10° 
e) 20 40 60 80 100 
TIME (ms) 


Figure 2.7 (a) Soot volume fraction, (b) particle diameter, and (c) number 


concentration along the line exhibiting the maximum volume fraction value. 


Only the curves labeled 3.85 cc/s are related to the ethene diffusion flame 


under investigation. Reproduced from ref. 86. 
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(c) type I probe © 


Figure 3.1 Sketch of thermophoretic probe devices; (a) circular grid, (b) bulk 


specimen carrier, (c)-(f) modified probes. 
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Figure 3.2 Sketch of thermophoretic probe sampling of a coannular diffusion 
flame. The flame coordinate system and the radial position r, of maximum 


soot volume fraction are also shown. 
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Figure 3.3 (A) Numerically calculated temperature contours, and (B) velocity 
contours of the flow field around the probe. The thermal boundary layer 6, 
and the trajectory of a soot particle ultimately captured halfway along the 
probe are also shown. (C) Schematic of a particle P traversing the distance 
between two neighboring isotherms. 
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Figure 3.5 Typical oscilloscope record of the probe trajectory. Transit time 


T,~16 ms, exposure (residence) time T,=36 ms. 


=162= 


RN SN NS NSP NS NGO NIIP SOREN Oe NGS 


Se VILLI ID FLT GG ED SUS ET TET EF ET (TAT UD EFT (TET ED OF LT 3 Ber 


Figure 3.6 Cross sectional drawing of the coannular burner. 
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Figure 3.7 TEM micrograph of soot particles deposited on a carbon substrate 
Supported by a copper grid positioned at z=30 mm. The grid bar is shown as 
the dark boundary. Grid exposure time T.=30 ms. 
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Figure 3.8 TEM micrographs of soot particles collected from the annular region 
of the ethene diffusion flame at various heights z, using carbon coated grids. 
For all cases T,=30 ms except (A) where T,=120 ms. (A) z=10 mm, (B) z=20 
mm, (C) z=30 mm, (D) z=40 mm, (E) z=50 mm, (F) z=80 mm. 


=L00— 


‘YIMOIS O[OTJIVd JOOS JO SWISTUBYOIP, GE 91INBI-Y 


ae 


yMoly aepns 


YyIMoly adepNsS uoHe1WO[S3Y 


-166- 


Figure 3.10 TEM micrographs of soot particles collected from the annular 
region of the ethene diffusion flame.at a height z=20 mm, using carbon coated 


grids. (A) T,=30 ms, (B) T,=50 ms, (C) T,=70 ms, (D) T,=90 ms. 
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Figure 3.11 High magnification 
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Figure 3.12 TEM micrographs of soot particles collected from the methane 
coannular diffusion flame using the type II probe. (A) z=25 mm, r=4.03 mmer,, 
T,~64 ms, (B) z=25 mm, r=4.77 mm>r,, T,~64 ms, (C) z=15 mm, r=4.43 mmer,, 


T,=72 ms, (D) z=15 mm, r=5.44 mmp>r,, T,=72 ms. 
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Figure 3.13 Primary particle size as function of height, near f=1 fore 
ethene diffusion flame. Each size measurement is an average over about 100 


primary particles. 
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Figure 3.14 TEM micrographs of soot particles collected on a type I probe 
from an axial location z=30 mm of the ethene diffusion flame, and at three 
different radial positions (A) r=1.28 mm, (B) r=2.54 mm-er,, (C) r=2.89 mm. 
Probe exposure time T,=36 ms. 
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Figure 3.15 TEM micrographs of soot particles collected on a type III probe 
positioned at z=10 mm in the ethene diffusion flame.Probe exposure time T,=48 
ms. (A) r=1.54 mm, (B) r=3.96 mmeer.C) r=4.59 mm. 
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Figure 3.16 TEM micrographs of soot particles collected from the fuel side of 
r(r<r.), at z=30 mm of the ethene diffusion flame. Type IV probe exposure 


time T,=30 ms. (A) r=0.64 mm, (B) r=1.30 mm. 
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Figure 3.17 TEM micrographs of soot particles collected from the outside 
region of r(r>r.), at z=30 mm of the ethene diffusion flame. Type IV probe 
exposure time T,=36 ms. (A) r=2.59 mm-sr,, (B) r=3.13 mm, (C) r=3.35 mm. 
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Figure 3.18 TEM micrographs of soot particles collected on a type IV probe 
positioned at z=15 mm in the ethene diffusion flame. Probe exposure time 
T,.=36 ms. (A) r=0.64 mm, (B) r=2.60 mm; (C) r=3.44 mm. 
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Figure 3.19 TEM micrographs of soot particles collected on a type I probe, 
positioned at z=l15 mm in the ethene diffusion flame. Probe exposure time 
T,=36 ms. (A) r=3.05 mm-r,, (B) r=3.73 mm, (C) r=3.05 mm~r,,(D) digitized 
picture of the particles shown in C. : 
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Figure 3.20 Comparison of the actual with the equivalent log-normal 
distributions at the two radial positions analyzed at z=15 mm. (A) r=3.05 mm, 
(B) r=3.73 mm. 


-177- 


25 


Computed Distr. Log-normal Distr. 


Ae iearphomyns & 


20 


15 
> 
: 
o > 
z 
o : 
® 3 
8 aa 
j : 
t : 
5 i 
acaba Y) 
5 0 wath: 
: a: 
g : mgs 
E sag? 
aa: 
; as 4: BA: 
15.0 35.0 55.0 75.0 95.0 
Diameter D (nm) 
150 
@ Log-normal Distr. 
4 Computed Distr. 
125 
J 
‘l 
100 i a: 
VY B 
4: 
> i B 
2 cae 
= 75 0 
a ‘4 
5 Matt 
ce ala 
. a ‘ 
aa: Bas 
50 (H ai - 
= Fy 3 - 
V a: ing: 
é it aracvedy 
ANAL a 
‘02: 0G: 8A: 8a: 
AYA VHA YE 
} 04 nA: 
AC 


15.0 35.0 55.0 75.0 95.0 
Diameter D (nm) 


Figure 3.21 Comparison of the distributions obtained from the numerical 


experiment described in the text. (A) 90 randomly distributed diameters, (B) 
500 randomly distributed diameters. 
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Figure 4.1 Schematic of a bimodal PVDF. The various processes involved in 


the formulation of the bimodal integral solution are shown. 
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Figure 4.4 Comparison of f, profiles obtained by the three methods, for the 


case of free molecular coagulation with a Gaussian source term. 
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Figure 4.5 Mass concentration histograms obtained by the MAEROS solution for 
the case of free molecular coagulation with a Gaussian source term. (A) t=0.02 


sec, (B) t=0.04see, (C)et=006tsec. 
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Figure 4.6 Particle volume distribution obtained by the two integral solutions, 
at an imstant where the source term achieves a maximum (t_,=0.03 sec). tie 


geometric mean volume of the source particles is also shown. 
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Figure 4.7 Time evolution of the PVYPF of the growth mode for the case of 


free molecular coagulation with a Gaussian source term. 
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Figure 4.9 Time variation of the geometric standard deviation of the growth 
mode for the case of free molecular coagulation with a Gaussian source term. 


The asymptotic value corresponding to the self-preserving form is also shown. 
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Figure 4.11 (A) Particle diameter and (B) number concentration versus 
residence time, as calculated by the bimodal solution in the numerical 
Simulation of the soot aerosol. The values determined by the optical and 
thermophoretic techniques are also shown. 
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Figure 4.12 (A) Time variation of fy and (B) Gaussian source pulse used in 


the numerical simulation of the soot aerosol. 
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Figure 4.14 Time variation of (A) soot surface area and (B) specific growth 


rate used in the numerical simulation of the soot aerosol. 


=193= 


fom 
‘s : —— With Surface Growth 
eC; lar seceeeees W/o Surface Growth 
Pe 
S 
2 12 
See 
= 
Ws 
2 T 
<= 10 
Oo 
> 
= 
10 
4 10 
= z=S5mm = z=10mm z=20mm z=30mm 
= 9 
10 : 
Q.00 0.01 0.02 OROR: 0.04 OOS 0.06 
(sec 
200.0 ( ) 
—— 
S 
Sf 5900 —— With Surface Growth 
a reeeeee W/o Surface Growth 
< 
au 
o 
100.0 
5 
Oo (B) 
Oo 
—b 9950.0 
= 
@ 
he 
0.0 


Figure 4.15 Time variation of (A) particle number concentration, and (B) 


particle diameter with and without the surface growth reaction term. 


=494> 


NBS-114A (REV. 2-8C) 
U.S. OEP T. OF COMM. 1. PUBLICATION OR 2. Performing Organ. Report NoJj 3. Publication Date 


BIBLIOGRAPHIC DATA REPORT NO. 


SHEET (See instructions) NBS/GCR-87/532 December 1987 
TITLE AND SUBTITLE 


Thermophoretic Sampling and Soot Aerosol Dynamics 
of an Ethene Diffusion Flame 


5 


« AUTHOR(S) 
Constantine M. Megaridis 


6. PERFORMING ORGANIZATION (If joint or other than NBS, see instructions) 


7. Contract/Grant No. 


8. Type of Report & Period Covered 
May 1986 


Brown University 
Providence, RI 


SPONSORING ORGANIZATION NAME ANDO COMPLETE ADORESS (Street, City, State, ZIP) 


National Bureau of Standards 
Department of Commerce 
Gaithersburg, MD 20899 


10. SUPPLEMENTARY NOTES 


[__] Document describes a computer program; SF-185, FIPS Software Summary, is attached. 


11. ABSTRACT (A 200-word or less factual summary of most significant information, If document includes a significant 
bibliography or literature survey, mention it here) 


A detailed investigation of the sooting structure of an overventilated laminar 
axisymmetric ethene/air diffusion flame under atmospheric pressure is presented. A 
thermophoretic sampling technique using a variety of probes is employed to obtain data 
on the morphological character of soot particles collected from various locations of 
the flame. The morphological features of the particles examined by Transmission Electra 
Microscopy provide not only qualitative information on particle agglomeration, surface 
growth, and oxidation but also quantitative data on primary particle size as a function 
of flame coordinates. An increasing primary particle size with height above the burner 
mouth proved the dominance of surface growth reactions on the morphology of particles 
transported in the lower and intermediate portions of the flame. The primary particle 
Size was found to decrease as the particles move to the upper regions of the flame 
probably because of oxidation. A theoretical model examining the aerosol dynamics 

when particle formation, coagulation, and surface growth are present is also discussed. 
The combined observations of the thermophoretic method and optical measurements by 
others led to the conclusion that the soot particles are formed lecally low in the 
flame in regions of high temperature near the reaction zone. 


12. KEY WORDS (Six to twelve entries; alphabetical order; capitalize only proper names; and separate key words by semicolons) 


aerosols; diffusion flames; morphology; particle size; soot; soot formation 


[13. AVAILABILITY 


14. NO, OF 
PRINTED PAGES 
K_] Unlimited 


[_] For Official Distribution, Do Not Release to NTIS 


[_] Order From Superintendent of Documents, U.S. Government Printing Office, Washington, 0.C. 
20402. 


210 


15. Price 


i. Order From National Technical Information Service (NTIS), Springfield, VA. 22161 


$24.95 


USCOMM=sOC 6043-P 80 


ee ee - Ot ES ie pyr s BRP 


~- ee ee ae 
ais nae GUT o& Lab teomeR vagy gainer? 


a ’ = Ae ua 
— ~ 2°. — a yee A 
c 3 
“ 
eclmanvd [ororsh 2002 bns, Ti sat ois szodqare 
Pn oa ~ 
Sonal? eens He Breil" 
“ 4 ad 7 oie § a ; 
w — - te — —_— eS Se. ee A Oe AA ATS AEE LE LLL AML A a a ay 
elibiuspeM .M Shleaaee 
- oe - wee eee nN NR A A I CT AA tee oe de 
sf ta) Set « 7B wv av 545 6 PID: bi HQITAS (RACRE: 
- jx NFS8 
ta 
eer. wre wits i 3° 3% vale BMAV MOTIASIM 
byaebrie. io Os uf 
: " 7 3 ASR 
M Seed azatty so 
- Pian len a Sry in ne 
~oR Fr? 
1 25 . = 7 a" 
oN : 7” asst 
A, Lc mye arverne as 
: brit te a wu? set Bei sg % : 
' P he VE afer 12), LPL aie , As = 
7 i wontvooe eds 26 BOE Sugize vit 
rt et tepttib the\scoedsey 
gotey suptarios rt lqmpea Ske 
y Joos 30 143 8 taoreolee 
odal 29 %1 : LG reolodagrom ? 
; 9 etileup vino 368 Shave 
i ty 2 ut, (GolLjgaas 
‘ / 
4 c ss ¥ olny 
r ors } ewol 
09 Ssvo & 7 n3 ‘ 2 © 958D 
& BE T&S 7 7 ZO 9 c 
>» 1s. (£3wowts SOR tue D5 wolisivgnes ,folzemtali @ 
T! Leos hon bode 229 -ormd? ef9 io agotievaesdO 
4 e443 we af2 secs sobsulones ef] Cae 
a3 ’ ‘ teen ssiosvegme2 dgid to enaky 
~— _ ree epcihcatiengtetiageerivsalhte tie ——— OFEE LETHE. ete Sa le ramets 
> bs ; — heen aq igp.es estas Aetre lovwedoriqie 199! sine evipw) oF nig) 
A yhbe” [Weta [Sara sie irre” Cypotorigrca Te orutcl ered Bir he 
va 
we , 3 omer’ 0.05 
m — a cael a I ee eee ++ EG Alten ent ATE A a a 
5 ? 
é Sn Sl J 


Z - ooh 
: . | ication of (US) particttTnenetrer foetedh opie g dee 
on SD Rotperrwnl ,esthO greed inemawrod .2,U ,2ioemugeO Vo see brake nie gad 5 
a | Lee iz w Syut tive guriszce growth reaction term. 
\ ' 
Laltt ,AY loteleniva® (2IT 4) eotvee® scbamotni ispineoe T 


# ey 

> 

pi -_ 

<p tec A Aa Og A GN COO A A ALLL ILE 


Og s-+a07 DG perOaGy 


L te + - re 4 ] 
op : a ey 
: iD ee 

7 = cy 


: ~ > ee 1" 
< ea ts A 
B ais een, 
7 


a 
8 


